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ABS TRACT 

Recent  progress  in  the  fabrication  of  low  attenuation  optical 
fibers  raises  the  attractive  possibility  of  optical  communications  via 
these  waveguides.  Integration  of  different  optical  functions  on  one 
substrate  in  the  form  of  optical  circuits  will  form  ideal  terminals 
for  such  a  communication  network.  We  shall  discuss  the  theory, 
fabrication  and  performance  of  a  few  "building  stones"  of  such  optical 
circuits.  Optical  directional  couplers  capable  of  coupling  light  from 
one  channel  to  another  will  be  described  and  their  use  as  switches  will 
be  discussed.  Distributed  feedback  for  integrated  lasers  will  be 
analyzed.  Laser  action  in  a  GaAs  waveguide,  in  which  distributed 
feedback  was  supplied  by  corrugations  fabricated  on  the  surface,  will 
be  described.  Longitudinal  and  transverse  mode  control  of  such  a  laser 


will  also  be  discussed. 
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I.  INTRODUCTION 


One  of  the  first  exciting  applications  envisioned  for  the 
laser  a  decade  ago  was  optical  communication,.  The  enormous  information 
carrying  capability  of  laser  beams  can  be  well  understood  by  noting 
that  at  a  wavelength  of  lp  the  frequency  is  3 •10'’  GHz. 

This  application  has  not  materialized  yet,  largely  because  at 
the  high  optical  frequency  the  laser  beam  is  greatly  attenuated  in 
passage  through  a  cloudy  and  rainy  atmosphere.  Optical  waveguides  in 
the  form  of  small  diameter  (2-3  thousandths  of  an  inch)  glass  fibers 
have  been  available  for  quite  some  time.  However,  the  attenuation  of 
the  light  propagating  in  those  flexible  fibers  was  so  great  that  their 
use  as  "optical  cables"  was  ruled  out. 

Recently  low  loss  optical  fibers  have  been  fabricated.  The 
fibers  can  be  used  for  short  distance  conmunication  systems  like  the 
ones  between  a  computer  and  its  peripheral  equipment,  or  for  medium 
and  long  distance  communication  systems  with  a  repeater  every  several 
miles.  They  combine  small  size  and  light  weight  together  with  large 
bandwidths  (several  G  bit/sec  for  a  single  mode  fiber  and  ceveral 
M  bit/sec  for  a  multimode  fiber).  Interference  and  pick-up  problems 
are  easily  eliminated  by  coating  the  fiber  with  opaque  material,  and 
ground  loop  problems  do  not  exist  in  fiber  optics  systems. 

A  review  of  the  principles  of  operation  of  these  fibers,  their 
properties  such  as  attenuation  and  bandwidth,  and  their  input  and 


output  devices  is  given  by  S.  Scmekh  and  A.  Yariv^.  The  enormous 
advances  in  the  fabrication  of  low  loss  fibers  can  be  appreciated  by 
considering  the  spectral  attenuation  curve  (Fig.  1.1)  for  a  multimode 
fiber  fabricated  by  researchers  at  Corning  Glass  Works^.  There  are 
two  Important  low  loss  attenuation  spectral  regions  with  loss  as  low  as 
4  db/km  around  0.85u  (aluminum  gallium  arsenide  laser  and  light 
emitters)  and  1.06p  (neodynium  YAG  laser). 

The  low  attenuation  in  these  fibers  together  with  their  large 
bandwidth,  revive  again  the  hopes  that  optical  communications  will  play 
an  important  role  in  the  telephone  network  and  in  other  types  of 
communication  networks  such  as  the  ones  on  aircrafts  and  ships,  and  the 
ones  between  a  computer  and  its  peripherals.  Another  potential  appli¬ 
cation  for  fibers  may  very  well  be  cable  television. 

These  applications  of  the  optical  fibers  call  also  for  small, 
rugged  and  inexpensive  terminals  in  which  light  will  be  generated  and 
detected,  and  in  which  information  will  be  imposed  and  extracted  at 
high  rates.  It  is  also  reasonable  to  believe  that  features  such  as 
switching  and  filtering  of  optical  signals  will  have  to  be  carried  out 
at  these  terminals.  These  requirements  have  stimulated  the  work  in 
the  field  known  as  "Integrated  Optics." 


Spectral  attenuation  of  a  low  loss  multimode  optical  fiber 
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1.2  Integrated  Optics 

A  typical  conventional  "optical  circuit"  ia  described  in  Fig.  2. 

It  consists  of  a  light  source,  a  modulator,  an  electronic  driver  for 
the  modulator,  polarizers  and  a  lens .  It  is  bulky,  heavy,  very  often 
unreliable  and  expensive.  In  addition  the  interface  between  the  solid 
state  electronics  and  optics  is  complicated  and  usually  limits  the 
performance.  The  desire  to  overcome  the  same  kind  of  problems  as  the 
ones  faced  by  electronic  circuits  a  decade  ago  has  led  to  the  development 
of  a  new  research  field  called  "integrated  optics"^”*^.  The  idea  is 
to  incorporate  on  a  small  substrate  all  the  optical  components  required 
to  make  up  the  various  optical  circuits .  Lasers,  modulators  and  detectors, 
along  with  waveguides  {the  "optical  wires")  and  the  analog  components 
to  lenses,  prisms,  mirrors,  polarizers,  etc.  are  to  be  fabricated  on  a 
single  chip  yielding  small,  rugged,  reliable  and  inexpensive  optical 
circuits. 

Let  us  examine,  for  example,  some  of  the  things  that  can  be 
fabricated  on  a  chip  of  GaAs.  This  material  and  its  alloy  A1  Ga.  As 

X  X-X 

(k) 

are  highly  suitable  for  integrated  optics  purposes'  .  Waveguides  can 
be  easily  constructed  in  GaAs;  injection  lasers  are  made  of  this  material 
and  a  p-n  diode  when  reverse  biased  is  capable  of  detecting  light. 

The  electrooptic  effect  (change  in  the  index  n  with  the  applied 
electric  field)  in  GaAs  is  among  the  largest,  thus  allowing  the  use  of 
efficient  modulators.  Finally,  the  possibility  of  incorporating  fast 
electronic  circuits  together  with  optical  circuits  on  the  same  chip  would 
ease  the  interfacing  problem. 
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Since  the  field  of  integrated  optics  is  in  its  first  stages, 
most  of  the  work  in  the  field  is  concentrated  on  providing  the 
building  stones  of  the  optical  circuits,  ftiese  building  blocks  include 
components  and  devices  such  as  guides,  directional  couplers 
modulators^,  detectors^,  surface  lasers  and  cavities ^1®)  mirrors, 
polarizers and  so  on. 

Chapters  II  and  III  deal  with  the  theory,  fabrication  and 
performance  of  some  of  these  components  and  devices.  Figure  1.3 
describes  two  types  of  dielectric  channel  waveguides  to  be  discussed. 
The  one  on  the  left  is  an  embedded  guide  while  the  other  is  a  ridged 
guide.  When  two  guides  are  closely  spaced  so  that  the  mode  profiles 
overlap,  optical  coupling  or  tunneling  occurs.  Figure  1.4  describes 
such  a  case.  If  the  interaction  length  is  appropriate,  light  entering 
the  device  in  one  channel  will  completely  tunnel  into  the  other  channel 
and  will  emerge  from  there.  This  device  iB  known  as  the  directional 
coupler,  and  we  shall  describe  the  first  experimental  demonstration  of 
the  effect  in  embedded  and  ridged  channels.  Controlling  the  tunneling 
or  the  coupling  between  the  two  channels  by  applying  electric  field 
turns  the  device  into  a  switch,  capable  of  steering  light  pulses  from 
one  channel  to  the  other  on  demand.  This  may  find  an  important  use  in 
optical  switching  and  multiplexing  systems.  We  shall  discuss  this 
possibility  in  detail  and  point  out  a  way  of  performing  it  with  lower 
electric  fields. 

In  our  treatment  of  the  theory  we  shall  attempt  to  avoid 
elaborate  solutions  of  Maxwell  equations,  but  rather  resort  to 


approximations  that  will  enhance  the  physical  understanding . 

In  Chapter  III  we  shall  concern  ourselves  with  the  linear 
applications  of  periodic  structures  in  thin  films.  Our  main  interest 
will  be  in  the  formation  of  a  suitable  cavity  for  an  integrated 
laser.  The  conventional  GaAs  laser  is  cleaved  on  both  ends.  The 
cleaving  forms  two  distinct  partial  mirrors  which  have  sufficient 
reflection  to  sustain  the  oscillation.  This  however  is  not  compatible 
with  the  concept  of  integrated  optical  circuits,  and  alternatives  will 
be  examined.  We  shall  describe  a  novel  alternative  which  uses  Bragg 
reflection  from  a  mechanical  corrugation  on  the  surface  of  the  wave- 
guide  as  shown  in  Fig.  1.5.  When  the  corrugation  period  A  satisfies: 


where  X  is  the  wavelength  of  the  guided  propagating  mode  6^  ,  the 
Bragg  condition  is  fulfilled  and  backward  reflection  occurs.  The 
intensities  of  the  forward  and  backward  propagating  modes  along  a 
distributed  mirror  of  length  L  are  shown  in  Fig.  1.5. 

For  GaAs  (Aq  =  0.83p  ,  n  =  3.6)  a  corrugation  with  a  very 
small  period  of  A  -  0.115C4  is  necessary  and  a  method  was  devised  to 
fabricate  it.  If  the  corrugation  extends  along  the  whole  length  of 
an  amplifying  medium  a  distributed  feedback  laser may  result.  The 
frequency  of  oscillation  of  this  laser  is  related  to  the  corrugation 
period  as  in  the  last  equation  and  can  be  varied  by  changing  A  . 

This  effect  is  potentially  useful  for  a  frequency  multiplexed  communi¬ 


cation  network. 


guide 


substrate 


Fig.  1.5  Bragg  reflection  from  a  distributed  mirror,  made  of 
corrugation  on  the  surface  of  the  waveguide.  The 
input  beam  decays  as  a  function  of  z  ,  giving  rise 


to  the  reflected  beam. 


We  shall  describe 
feedback  laser  in  GaAs^^ 
mode  control. 


the  first  observation  of  a  distributed 
and  elaborate  on  longitudinal  and  transverse 


1.4  Fabrication  of  Integrated  Optics  Components 

We  shall  devote  a  whole  chapter  (Chapter  IV)  zo  the  fabrication 
of  the  devices  we  have  discussed  earlier,  because,  in  our  opinion,  the 
fabrication  problems  will  determine  the  rate  at  which  the  new  field  of 
integrated  optics  will  develop.  The  dimensions  of  the  channel  guides 
are,  as  in  their  microwave  counterparts,  of  the  order  of  the  wavelength. 
This  means  that  structures  with  dimensions  of  a  few  microns  have  to  be 
fabricated  with  high  quality  edge  smoothness.  As  another  example  we 
recall  that  for  a  distributed  feedback  laser  in  GaAs  a  corrugation  with 
a  period  of  0.115CU  has  to  be  fabricated  on  the  surface  of  GaAs. 

In  collaboration  with  Hughes  Research  Laboratories  (Malibu,  CA), 
we  have  employed  ion  beam  milling  and  ion  implantation  in  the  fabrica¬ 
tion  of  the  different  devices.  We  have  also  devised  new  techniques, 

(12) 

such  as  holographic  exposure  of  photoresist  combined  with  ion  milling'  ' , 
for  the  fabrication  of  the  small  period  corrugation. 
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II.  OPTICAL  DIRECTIONAL  COUPIERS 

II , 1  Introduction 

The  optical  directional  coupler,  analogous  to  the  microwave 
element^ ^  of  the  3ame  name,  consists  of  parallel  channel  optical 
waveguides  sufficiently  closely  spaced  that  energy  is  transferred  from 
one  to  another.  For  this  coupling  to  take  place  cumulatively  over  a 
substantial  length,  the  light  must  propagate  with  the  same  phase 
velocity  in  each  channel.  The  fraction  of  the  power  coupled  per  unit 
length  is  determined  by  the  overlap  of  the  modes  in  the  separate  channels. 
Thus,  it  depends  on  the  guides'  separation,  the  mode  penetration  into 
the  substrate,  and  the  interaction  length.  Figure  II. 1  is  an  example 
of  a  multichannel  directional  coupler.  It  shows  a  diagram  of  a  large 
number  of  coupled  channel  waveguides  (produced  by  proton  implantation) 
and  typic«l  intensity  profiles  of  the  guided  lignt.  The  incident  light 
is  focused  into  a  single  channel  at  z  =  0  ,  but  is  coupled  into  the 
adjacent  guides  as  it  propagates. 

A  theoretical  analysis  of  the  coupling  between  dielectric 

(2) 

waveguides  usually  involves  a  formal  solution  of  Maxwell's  equations'  . 

We  shall  attempt  here  to  derive  an  expression  for  the  coupling  coefficient 
and  its  sign  from  an  intuitive  physical  point  of  view  without  resorting 
to  an  elaborate  mathematical  solution.  We  shall  also  derive  the  power 
distribution  as  a  function  of  propagation  distance  for  two  coupled 
guides  and  an  infinite  number  of  coupled  guides. 

Experimental  results  of  a  dual  channel  directional  coupler 
and  a  multichannel  coupler  will  be  given.  These  channels  were  embedded 
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at  the  surface  of  GaAs  by  proton  implantation.  A  different  type  of 
directional  coupler  suitable  for  ridged  waveguides  is  also  described. 
We  conclude  with  an  account  of  the  different  applications  of  these 
devices  in  an  integrated  optics  circuit,  such  os  couplers,  modulators, 
multiplexers  and  polarizers. 

11.2  Coupled  Mode  Formalism 

Let  us  consider  the  two  waveguides  shown  in  Fig.  11.2a.  Hie 
waveguides  can  have  an  arbitrary  cross  section;  however,  for  the  sake 
of  simplicity,  we  consider  only  guides  with  rectangular  cross  sections. 
The  electric  field  of  a  propagating  mode  in  a  guide  is  described  by 

-t£3l«9 

the  complex  amplitude  A(z)  which  includes  the  phase  term  e  p  "  ,  and 
by  the  mode  C(x,y)  solved  for  in  the  absence  of  the  other  guide. 

We  thus  have: 


E(x,y,z)  =  A(z)  6(x,  (II. 2.1) 

The  mode  profiles  e(x,y)  are  always  normalized  to  carry  one  unit  of 
power;  the  power  carried  by  E1(x,y,z)  (where  the  subindex  1  denotes  the 
guide  number),  at  a  given  z  is 

Px(z)  =  jA1(z)|2  =  Ax(z)  A*(z)  (II. 2. 2) 

(1) 

Employing  the  coupled  mode  formalism'  ,  we  write  the  general  coupled 
mode  equations  for  the  amplitudes  of  the  two  modes 


dA0(z) 

dz 


"i(Wz*  '  *oiAi(8) 


( 11.2. 3) 


n -ft  tySJW  '  *f 
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<JA.(z) 

-J—  =  -iPjA^z)  +  K^z)  (11.2.4) 

PQ  and  P1  are  the  propagation  constants  of  the  modes  in  the  two 
guides  (an  imaginary  part  in  p  describes  loss  or  gain  that  may  exist 
in  the  guides).  K01  and  represent  the  coupling  between  the  two 

modes,  and  a  great  deal  of  the  theory  in  this  chapter  will  be  devoted 
to  their  derivation.  Before  going  into  that,  let  us  consider  a 
theoretical  and  an  experimental  example  of  the  coupled  mode  behavior. 

11.3  Dual  Channel  Directional  Coupler- Theory  and  Experiment 

Consider  the  guides  shown  in  Figure  II. 2a.  We  assume  that 
the  guides  sire  identical  and  that  both  have  an  exponential  loss  cansteuit 
of  a  .  Thus  the  propagation  constant  of  the  guide  is  given  by 

P  =  Pr  -  i  f  (II. 3.1) 

As  we  shall  see  later  the  coupled  mode  equation  for  this  case  can  be 
written  as: 


dA_(z) 

—2 -  =  -ip  A  (z)  -  iKA.(z) 

dz  u 


(II. 3. 2) 


dA.(z) 

— = -  =  -iP  A.(z)  -  iKA  (z) 

dz  1  0 


(II. 3. 3) 


where  K  ,  the  coupling  coefficient,  is  a  real  quantity.  We  also 
assume  that  at  z  =  0  light  i3  coupled  into  guide  0  so  that  the 


Ac,  (xyz)  *  Ac  *  cQ  (n|  -  n|) 

Ac, (xyz)  *  0  elsewhere  _ w _ n3 

pig.  II. 2  (a)  Sketch  of  two  adjacent  rectangular  waveguides  and  their 

propagating  inodes. 

(b)  Definition  of  the  dielectric  discontinuity  Ae,(xyz). 
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boundary  conditions  for  the  problem  are  given  by 


A0(0)  =  1  A1(0)  =  0  (11.3.4) 

The  solution  is  described  by: 

AQ{z)  =  cos  Kz  ei3z  (II. 3. 5) 

A1(z)  =  -i  sin  Kz  ei3z  (11.3.6) 

The  power  flow  in  the  guides  therefore  is 

PQ(z)  =  AQ(z)  A*(z)  cos2Kz  e-02  (II. 3.7) 

P^z)  =  A1(z)  A*(z)  =  sin2Kz  e-0®  (11.3.8) 


From  II. 3. 7  and  11.3.8  we  learn  that  the  power  will  oscillate  back  and 
forth  between  the  two  guides  as  shown  in  Fig.  II. 3.  Because  of  tne 
losses  present  in  the  guide,  the  power  will  also  decay  exponentially  as 
a  function  of  the  propagation  distance.  The  length  necessary  for  a 
cample te  transfer  of  power  from  one  guide  to  the  other  is  given  by 

L  =  (II. 3. 9) 

Such  a  directional  coupler  that  allows  the  transfer  of  light 
from  one  channel  to  suiother  is  one  of  the  building  blocks  of  optical 
integrated  circuits.  We  shall  describe  in  what  follows  the  first 
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operation  of  such  a  device  in  the  optical  regime. 

The  dual  channel  directional  coupler  was  fabricated  by  proton 
implantation  in  GaAs in  a  fashion  described  in  detail  in  Chapter  IV. 
The  channel  guides  were  formed  by  implanting  the  protons  through  a 
gold  mask  with  the  appropriate  slit  openings  in  it.  Bombardment  causes 
defect  centers  in  the  material,  which  trap  the  free  carriers.  This 
compensation  of  the  free  carriers  eliminates  their  negative  plasma 
contribution  to  tV.  index  of  refraction,  and  thus  increases  the  index 
for  the  samples  used  here  (Nsuj,strate  ~  2  x  10  )  by  An  .005  . 

The  cross  section  of  the  gold  mask  on  top  of  the  GaAs  substrate  is  shown 
in  Fig.  II. Ua.  This  way  the  guides  embedded  at  the  surface  have  a 
cross  section  of  about  3x3 u  ,  an  index  discontinuity  of  about 
An  =“■  .005  ,  and  a  separation  of  3n  between  the  guides. 

Fig.  II. 5  describes  three  pair  of  such  directional  couplers. 

On  the  left  is  the  input  plane  of  the  device  and  the  black  dot  indicates 
the  point  at  which  the  input  light  was  focused.  On  the  right  is  a 
series  of  scans  of  the  light  distribution  at  the  output  plane,  taken  by 
the  set-up  described  in  Ref.  (3).  Let  us  confine  our  attention  to 
the  directional  coupler  at  the  center.  The  second  row  from  the  top 
shows  that  when  the  input  light  is  coupled  into  the  channel  on  the 
right,  it  emerges  at  the  output  plane  from  both  channels .  Similarly, 
the  second  row  from  the  bottom  indicates  that  the  light  also  emerges 
from  both  channels  when  it  is  coupled  into  the  left  channel.  Since  the 
intensity  of  the  emerging  light  from  both  channels  is  about  equal,  we 
can  write  for  this  sample, 


Fig.  II. 4  Scanning  electron  microscope  photograph  of  the  gold  mask  on  the 
GaAs  substrate  used  in  the  fabrication  of  the  (a)  dual 
directional  coupler,  (b)  multichannel  directional  coupler.  The 
remains  of  the  photoresist  can  be  seen  on  top  of  the  gold  stripes. 


mm  long 


Sketch  of  the  input  plane  of  a  1  mm  long  directional  coupler. 

On  the  right  are  photographs  of  guided  light  intensity  profiles 
at  the  output  for  various  input  conditions . 


-24- 


KL 


Jt 

T* 


t  =  1  mn 


K  *■  O.79  nm" 


By  increasing  the  length  of  the  sample  by  about  a  factor  of  two,  a 
complete  transfer  of  power  should  be  observed  because 


Iliis  is  shown  in  Fig.  11.6.  The  second  row  from  the  top  describes  a 
situation  in  which  light  coupled  into  the  left  channel  emerges  from 
the  one  on  the  left,  and  the  second  row  from  the  bottom  shows  the 
transfer  of  power  from  left  to  right. 

This  device  can  be  turned  into  a  switch  by  controlling  K  , 
the  coupling  coefficient,  with  the  application  of  an  electric  field 
through  the  electrooptic  effect.  This  exciting  possibility  will  be 
discussed  in  detail  later. 

Let  us  go  back  now  and  derive  theoretical  expressions  for  the 
value  of  the  coupling  coefficient. 

II. 4  Derivation  of  the  Coupling  Coefficient 

We  consider  again  the  coupled  mode  equations: 


dA  (z) 

-0 -  =  -iP0A0(z)  +  K^A^z)  (II.  2. 3) 

dz 


dA^z) 


dz 


=  -iPjA^z)  +  K10Aq(z) 


(II.H.4) 


recti onal  Coupler  -  2.1mm  long 


Fig.  II. 6  Sketch  of  the  input  plane  of  a  2.1  ram  long  directional  coupler. 

On  the  right  are  photographs  of  guided  light  intensity  profiles 
at  the  output  for  various  input  conditions. 
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We  also  recall  that  the  electric  field  in  a  guide  is  given  by 

E(x,y,z)  =  A(z)  e(xy)  (II.2.1) 

and  that  the  power  carried  by  E^y,*),  for  example,  is 

Pl(z)  *  lAx(z) 1 2  =  Ax(z)  A*(z)  (II. 2. 2) 


Next  we  obtain  the  rate  of  power  growth  in  mode  1  due  to  the  coupling 
from  mode  0.  (For  this  purpose  we  assume  that  p  is  real.)  To  do 
this  we  multiply  11.2.4  by  A^(z)  and  the  complex  conjugate  of  11.2.4 
by  A^(z)  .  This  gives 


*  dAj^z)  # 

Ai(z)  dz  “  *^iAi(z)  Ax(z)  +  kioA0(z)A1(z)  (n-4.1) 


Ax(z) 


dAx(z) 


*  * . 


dz 


=  13qA^(z)  A^(z)  +  K10Aq(z)A1(z)  (11.4.2) 


Summing  the  two  equations  thus  yields 


d[A.(z)  A*(z)] 

di -  =  2  RetA^z)  Kjq  AQ(z)  ]  (11.4.3) 

Using  II. 2. 2  gives 


dPl(z)  *  * 

“~dz  =  RefVz>  Kqi  V2)3  (11.4.4) 


-  - -  ■  ■ — ^-i^- 


- - 
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We  shall  now  obtain  a  different  expression  for  the  rate  of 
power  growth  in  guide  1  from  a  simple  minded  physical  point  of  view, 
in  the  following  manner.  Guide  0  in  the  absence  of  guide  1  supports 
a  mode  with  a  mode  profile  given  by  6Q(xy)  .  This  mode,  in  spite  of 
having  a  long  tail  in  the  x  and  y  direc-xons,  suffers  no  loss,  and 
propagates  with  a  constant  power.  When  gpHde  1  is  introduced,  tne  tail 
of  6Q(xy)  penetrates  into  guide  1  and  generates  additional  polarization 
because  of  the  higher  index  inside  the  guide.  This  perturbation  in  the 
polarization  in  turn  drives  the  mode  in  guide  1. 

The  net  power  per  unit  volume  expended  by  an  electric  field 
on  the  electric  polarization  is  given  by  the  veil  known  relation 


Power  _  p  5P 
Volume  “  dt 


(IX.U.5) 


where  E  is  the  field,  P  is  the  polarization,  and  the  horizontal  bar 
denotes  time  averaging. 

We  can  thus  express  the  power  generation,  due  to  the  perturba¬ 
tion  in  the  polarization  at  a  given  point  in  guide  1,  by  changing  the 
sign  on  the  right  of  equation  II. 4. 4  yielding 


dP,0(xyzt) 
=  -E1(xyzt)  •  - 


(II. 4. 6) 


where  P^q  is  the  driving  polarization  caused  by  the  mode  in  guide  0, 
This  perturbation  in  the  polarization  can  be  approximated  by  multiplying 
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the  electric  field  of  mode  0  by  the  disturbance  of  the  dielectric 
constant  caused  by  the  introduction  of  guide  1. 


=  EQ(xyz)  •  Ac^(xyz) 


(11.4.7) 


Over  the  cross  section  of  guide  1,  Ae^(xyz)  has  a  constant  value, 
independent  of  z  ,  of 


Ae  =  e2-e3  =  EQ(n2  -  n^) 


and  is  zero  elsewhere  as  shown  in  Fig.  II. 2b.  Using  II. 2. 3, 
II. U. 7  and  time  dependence  of  e*^  ,  II. 4. 6  becomes: 


—  -  -A1(z)ela5tei(xy)AE1(xy)iaA0(z)eia3te0(xy)  (II. 4, 5) 


Performing  the  time  averaging  yields 


r  =  ^(A^zJ^zyJC-laA^xyJC^xyjAQtz)]*}  (11.4.9) 


The  growth  of  the  total  power  P^  carried  by  the  mode  in  guide  1  can 
now  be  derived  by  simply  integrating  over  x  and  y  .  We  obtain 


dPx(z) 


iRe|A1(z)[icuJ’  ei(xy)Ae*(xy)£0#(xy)dxdy]A*(z)| 


(11,4.10) 
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This  can  also  be  written  as: 

dPi(^)  /  4  *  *  v 

- - -  ^RelA^zJC-ltoJ  ei(xy)A£l(xy)e0(xy)dxdy]  A0(z)J  (II.U.ll) 

A  comparison  of  II.U.U  and  II.U.ll  immediately  reveals  that 

Kio  =  nr  J  (n.u.12) 

Similarly  one  obtains 

Koi  =  I  ^(xy )A£Q(xy)61(xy )dxdy  (II.U.13) 

where  AcQ(xy)  as  opposed  to  Ae^xy)  ,  is  equal  to  Ac  over  the  cross 
section  of  guide  0  ,  and  is  zero  elsewhere. 

Fran  II.U.12  and  II.U.13  it  is  evident  that  as  long  as  A£ 
is  real  and  the  guides  are  similar, 

K10  -  -K0X  (II.U.1U) 

As  shown  in  Appendix  II  this  is  a  general  property  of  lossless  co- 
directional  coupling.  It  is  interesting  to  note  that  this  relation 
holds  also  for  lossy  propagation  in  which  0  is  complex  but  Ac  is 
real  (such  a  case  may  occur  when  the  absorption  in  the  guide  and  in  the 
substrate  are  the  same,  or  when  the  losses  are  caused  by  scattering). 

In  the  cases  we  consider  Ac  is  mainly  real.  In  general  a  conplox 
Ac  will  modify,  for  example,  the  solutions  described  in  Fig.  II. 3  for 
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two  modes,  by  causing  the  two  modes  to  carry  the  same  power  when  z 
is  sufficiently  large^ . 

Equations  11.4.12  and  11.4.13  are  identical  to  those  obtained 

(4) 

by  Marcuse  '  who  solved  the  Maxwell  equations  using  perturbation 
theory.  They  are  independent  of  the  guide  geometry  and  can  be  applied 
to  various  cross  section  guides  once  the  profile  of  the  propagating 
mode  is  known. 

The  problem  of  coupling  between  the  guides  can  be  approached 

from  a  different  point  of  view,  in  which  the  coupling  coefficient  is 

derived  by  solving  for  the  propagation  constants  of  the  two  eigenmodes 

of  the  two  guide  directional  coupler,  as  shown  in  Fig.  II. 7.  The  two 

modes  of  a  symmetric  coupler  are  a  symmetric  one  £  and  an  anti- 

s 

sym  trie  one  £  ,  with  propagation  constants  S  and  p  .  If  at 
z  =  0  light  is  coupled  into  the  0  guide  only,  the  two  modes  are 
excited  equally  in  such  a  way  that  they  add  on  the  0  guide  side  and 
cancel  over  guide  number  1  side.  Since  B  ^  p  the  two  modes  drift 
out  of  phase  as  they  propagate  along  the  z  direction.  At  a  distance  L 
which  is  given  by 

AfJL  =  (0s  -  0a)L  =  *  (11.4.15) 

the  two  modes  are  180C  out  of  phase,  cancel  over  guide  number  0  and 
add  over  guide  number  1,  as  shown  in  Fig.  II.7.  The  net  result  is  a 
complete  transfer  of  power  from  guide  0  to  guide  1.  From  a  comparison 
of  11.4.15  and  II. 3 .9  we  find  that  AP  is  related  to  the  coupling 
coefficient  by 


K 


(II. 4. 16) 


guide  0 

guide  1 

Fig.  II. 7  The  symmetric  and  antisymmetric  eigen  modes  of  a  dual 
directional  coupler.  Power  is  transferred  from  one 
channel  to  another  because  of  not  identical  propagation 
constants. 
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and  nQ  is  the  index  of  refraction  distribution  in  the  absence  of 
guide  1,  while  is  the  distribution  in  the  absence  of  guide  0.  For 
effective  exchange  of  power  between  the  two  guides  the  propagation 
constants  must  be  similar.  In  section  II. 10  we  shall  treat,  from  a 
coupled  modes  point  of  view,  a  case  where  they  are  not  similar.  Let  us 
therefore  confine  our  attention  now  to  a  case  where 


P0  =  pi  s  ^  (11.4.19) 

As  will  be  discussed  in  section  II. 6  we  can  assume  that  the  modes  have 

one  main  component  of  the  electric  field,  say  in  the  y  direction,  and 

therefore  a  mode  can  be  normalized  to  carry  one  unit  of  power  by 

imposing  (see  Appendix  I): 

00 

2^7  ||ex(xy)j2dxdy  =  1  i  =  0,1  (11.4.20) 

We  now  approximate  the  modes  of  the  two  guide  directional 
coupler  by  the  sum  and  difference  of  the  individual  modes. 


-  7-5 — T2TT  [aeo(iy)  4  ^l^)1 

(a  +  b  }b 


C  (xy)  =  — - — -  [bC  (xy)  -  a£  (xy)] 

(a2  +  b2)fc 


(11.4.21) 


To  ensure  the  orthononnality  of  these  two  modes  we  have  to  assume  a 
small  overlap  of  the  individual  modes. 

J|ei(xy)CQ(xy)jdxdy  «  J  |C0(xy) |2dxdy  (11.4.22) 

The  orthogonality  of  £  and  £  is  important  because  it  ensures  the 

S  cl 

conservation  of  energy,  once  a  combination  of  these  modes  is  excited 

and  propagates  in  the  z  direction. 

To  find  the  propagation  constants  of  the  above  modes,  (3  and 

s 

insert  them  into  the  appropriate  wave  equation. 

,  CD  2  p  fl2p 

V. 4  7  aiA  ■  p.es 


7 \i  +  2L  4-e  .  p2e 
t  a  2  II  a  &  a 
c 


(11.4.23) 


where  is  the  index  distribution  when  both  guides  exist  etnd  i3 

shown  in  Fig.  II. 2a. 

Let  us  consider  the  equation  for  £  .  With  the  aid  of 

s 

11.4.17  and  11.4.21  we  can  write  the  upper  equation  in  11.4.23  as 

“2^nil  *  n0^0  +  ""^nll  "  ni)^i  =  (^s  "  +  b^l^ 

c  c 


We  recall  that  (n^  -  n2)  and  (n2.,  -  n2)  are  actually  the  perturba¬ 
tions  in  the  dielectric  constant  we  have  discussed  earlier.  We  can  thus 


write 


(11.4.25) 
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nII  -  n0  *  V£o 
nII  *  n\  *  4e(/Eo 


where  Ae,  is  shown  for  example  in  Fig.  11.2b  and  is  the  free 

space  dielectric  constant.  Next  we  obtain  two  expressions  by  multiply 
*  * 

ing  11.4.24  by  t q  or  C1  and  integrating  over  x  and  y  .  Using 

11.4.25  and  c2  =  ~-  we  get 
M  0 

00  00  CD 


co^bje^Ae^dxdy  +  Aaje^Ae^dxdy  =  (02-p2  jJcjJ^o+bejJdxdy 

_  00  —00 

(11.4.26) 


-CD 


2 

u>  (ia 


>J4 


Ae^dxdy 


2  !  * 

co  M-oj  ^jAeQCjdxdy 


(Ps'p2)Jei(ae0+bei)dxdy 

(11.4.27) 


The  second  integrals  on  the  left  are  second  order  terms  because  As^  , 
for  example,  is  different  from  0  only  where  |  £Q  | 2  is  very  close  to 
zero.  Thus  using  11.4.22  we  can  write  11.4.26  and  11.4.27  as 


o>2 nbj £*A  e  jf^dxdy  =  a(02-p2)je*eodxdy  (11.4.28) 


u)2Maje*AEieodxdy  =  b02-62)j  e^dxdy  (11.4.29) 

.so  _co 


Taking  the  ratio  of  the  two^equations  we  get 

2  je^e^dxdy 


.00 

-w 


(11.4.30) 


Jo 


eiA'oeodx4,r 


In  a  similar  fashion,  starting  with  the  equation  for  in  II . t .23, 
we  find 


•  iiftt  tfa  ft  fa  Mt ' ( 
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V, 
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b2  jcjAe^dxdy 


(II.lt.31) 


a& 


J*l4'oPodx4lr 


From  11.4.30  and  II. 4.31  we  can  conclude  that  far  the  existence 
of  a  solution  In  the  form  ox  11.4.21  we  must  require 

00  CO 

Je^-^dxdy  =  j£*deoeodxdy 


(11.4.32) 


which  yields 


a  =  b 


(II.4.33) 

We  can  now  obtain  B  by  adding  11.4.28  and  II.4.2Q.  We  also 

S 

assume  that  B  -B  «  B  1111(1  with  the  help  of  11.4.20  we  get 

s  m  os 

=  B  tJeoA€oeidxdy+j ^i^Ci^Qdxdy ]  (11.4.34) 


In  the  same  fashion  we  obtain  for  B 


a 


=  ‘  ?JC^e0Cl<h^JClAelC0dx<,y]  (11.4.35) 


-00 


The  coupling  coefficient  can  now  be  obtained  from  II. 4. 16, 
11.4.34  and  11.4.35 

00  00 

K  “  ’  I  fJl'R04E0eldxdy+jei4EieodxdyI  (“•‘*•36) 

-OB  —OB 

The  value  of  K  is  in  agreement  with  the  expression  derived  earlier. 

It  is  reasonable  to  assume  that  II. 4. 36  yields  a  good  approximation 
for  K  also  in  a  case  where  the  guides  are  not  identical  and  II. 4.  32 
is  violated. 
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It  is  interesting  to  note  that  from  the  above  point  of  view  it 
in  possible  to  easily  obtain  an  upper  limit  for  the  value  of  the 
coupling  coefficient  K  .  Clearly,  this  value  reaches  its  upper  limit 
as  the  two  guides  cure  brought  closer  and  closer  together.  Eventually 
they  form  one  guide  with  twice  the  width.  The  symmetric  mode  turns 
into  the  lower  order  mode  of  the  wide  guide  and  the  antisymmetric  mode 
turns  into  its  higher  order  mode.  The  calculation  of  the  propagation 
constants  of  the  new  modes  is  straightforward  and  yields  the  upper 
limit  for  . 


II.  5  Coupling  between  Planar  Guides 

As  a  first  example  we  apply  our  result  to  the  planar  waveguides 
shown  in  Fig.  II. 8.  Ihe  mode  profile  of  a  single  planar  waveguide  can 
be  derived  analytically  (see  Appendix  I).  We  consider  a  TE  mode  which 
has  its  electric  field  in  the  y  direction,  carries  a  unit  power  per 
unit  length  in  the  y  direction,  and  is  given  byv  ' : 


where  kQ  is  the  free  space  wave  number  and 


hx  m(4kl  - p2)i 


•x  ■  ^ 


(n.5.2) 
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Fig.  II. 8  Sketch  of  two  adjacent  planar  guides  and  their 
propagating  modes. 
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In  addition,  tbe  eigenvalue  equation  imposes  the  relation  for  the  mth 
order  mode  (m+1  is  the  number  of  extrema) 


-1 

h  v  =  (m+l)jt  _  2  tan” 

x 


(II- 5. 5) 


A  similar  profile  describes  the  mode  in  guide  number  1,  the  only 
difference  being  that  the  center  of  the  mode  is  at  x  =  w+s  rather 
than  x  =  0. 

Knowing  the  mode  profiles  we  can  obtain  the  coupling  coefficient 
between  the  guides  with  the  help  of  equation  11.4.12.  We  note  first  that 
Ae^  for  this  case  is  given  by 


Ae^xy)  «  e2  -  e 


co^nl'nP  f0T 


s  +  ^xss  +  |! 

_00  S  y  S  _00 


(II. 5.4) 


AEi(xy)  =0  ,  elsewhere 


We  thus  obtain  the  following  expression  for  the  coupling  coefficient; 


(II. 5. 5) 


In  the  derivation  of  this  expression,  the  y  integration  is  performed 
over  one  unit  length  and  we  have  assumed  that 

jp  cos  h  ^  +  h  sin  h  »  |p  cosh  -•  -  h  sinh  ^|e  *** 
l4x  x  2  x  x  2'  ,rx  x  2  x  x  2 ' 

(II. 5. 6) 
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which  is  always  true  for  well  confined  modes. 

Equation  II. 5. 5  is  the  same  as  approximate  expressions  for  the 
coupling  coefficient  between  planar  waveguides  given  in  Refs.  4,5,6. 

In  Ref.  5  a  comparison  of  this  expression  to  the  exact  solution 
of  the  problem  is  made  and  a  good  agreement  is  shown. 

II. 6  Coupling  between  Channel  Guides 

A  rigorous  solution  for  the  mode  profile  at  a  rectangular 

(7) 

cross  section  channel  waveguide  requires  a  computer' 1 ' .  However,  once 
the  mode  profile  is  known,  11.4.12  can  be  applied  for  the  calculation 
of  the  coupling  coefficient.  Marcatili'  '  has  shown  that  it  is  possible 
to  introduce  a  drastic  simplification  in  the  derivation  of  the  mode 
profile  which  enables  one  to  get  a  simple  solution.  Consider  the 
waveguide  shown  in  Fig.  II .9.  Most  of  the  power  propagates  in  the  high 
index  region  and  the  fields  decay  exponentially  in  the  x  and  y  directions. 
Consequently,  for  a  reasonably  confined  mode  only  a  small  part  of  the 
power  propagates  in  the  four  shaded  areas.  This  means  that  in  matching 
the  boundary  conditions  the  field  along  the  edges  of  the  shaded  area 
can  be  ignored. 

The  solution  of  the  simplified  problem  yields  modes  that  are 

essentially  of  the  TEM  kind  and  can  be  grouped  in  two  families,  e* 

mn 

and  £?  .  The  first  family  has  most  of  the  electric  field  in  the  x 

mn 

direction  and  the  magnetic  field  in  the  y  direction,  while  the  second 
family  has  most  of  the  electric  field  in  the  y  direction  and  the  magnetic 
field  in  the  x.  Our  notation  is  such  that  ra+l  and  n+1  indicate 
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the  number  of  extrema  of  the  fields  in  the  x  and  y  directions.) 
Furthermore  the  solution  of  the  problem  indicates  that  mode  profiles  in 
the  x  and  y  directions  are  actually  given  by  two  equations,  each  of 
which  is  identical  to  the  eigenvalue  equations  one  obtains  for  a  one- 
dimensional  (planar)  waveguide.  For  example,  we  consider  an  $ 
mode  in  the  guide  shown  in  Fig.  II. 9.  The  eigenvalue  equation  for  the 
x  direction  is  identical  to  that  of  a  IE  mode  in  a  synmetric  planar 
guide  (see  Appendix  I).  It  is  given  by 

-1  k 

h  w  =  (m+l)rt  -  2  tan”  _JS  (II. 6.1) 

x  p 

*x 

where 

h2  +  =  (n2-n2)k2 

x  x  '  2  3'  0 


The  eigenvalue  equation  for  the  y  direction  is  identical  to 
that  of  a  TM  mode  in  an  asymmetric  planar  guide  (see  Appendix  I)  and  is 


given  by 


2  2 

,  n  h  1  nl  h 

h  t  =  (n+l)it  -  tan”  .  tan”  ? 

7  "2  Py  "3  -V 


(11.6.2) 


where 


.  2  2  .  2  2\.2  .2  2  /  2  2  \  .2  /  tt  \ 

hy  +  py  =  ^n2”n3)ko  ;  hy  +  V  =  (n2'n!)k0  (11-6.3) 


The  propagation  constant  along  the  z  direction  is  therefore  given  by 


(II. 6. 4) 


We  can  therefore  approximate  the  propagating  mode  by  regarding 
only  the  main  components  of  the  fields  and  using  the  mode  profiles  as 
given  by  II. 6. 2  and  II. 6. 3.  The  validity  of  this  approximation  can  be 
determined  by  comparing  the  value  of  the  expression 


according  to  the  figures  in  Ref.  2,  the  approximation  is  within  a  few 
percent  of  the  exact  value.  Equation  II. 6. 6  also  indicates  what  modes 
are  supported  by  the  guide,  because  only  modes  for  which  the  expression 
is  positive  are  above  cutoff. 

Figure  II. 9  shows  the  mode  profile  of  the  mode.  This 

mode  profile  after  the  appropriate  normalization  can  be  used  in 
conjunction  with  11.4.12  to  calculate  the  coupling  coefficient. 

We  shall  attempt  to  get  a  good  estimate  for  the  coupling 
coefficient  without  going  through  the  above  process.  We  note  that  the 
profile  of  the  mode  in  the  x  direction  is  the  same  as  that  of  the  planar 
guide  we  have  dealt  with  before.  The  major  difference  is  that  the  planar 
guide  and  its  propagating  modes  are  independent  of  y.  This  basically 


means  that  ot^(xy)  (which  is  defined  for  the  planar  guide  case  by 
II. 5.  *0  has  the  same  value  for  any  value  of  y  to  which  the  modes  of 
the  guides  extend.  IMs  would  also  be  the  case  in  a  channel  waveguide 


if  the  modes  were  very  well  confined  in  the  y  direction 
In  this  case  Ae^xy)  is  defined  as 


(i  ,i-  « i) 

Vt  py  / 


Af-^xy)  =  -  e  2  =  £Q(n^  -  n^)  for  s  +  |  •- 


x  -  s 


2 


■  t  y  0 


Aex(x y)  =  C 


elsewhere 


(11.6.7) 


Thus  if  the  mode  is  very  well  confined  in  the  y  direction  it  extends 
only  from  y  =  0  to  y  =  -t  which  is  the  region  at  which  Af(xy) 
has  a  constant  value.  For  this  case  II. 5. 5  can  be  used  to  express  the 
coupling  between  the  two  channels.  We  can  modify  II. 5. 5  to  include 
coupling  between  modes  that  are  not  very  well  confined,  but  satisfy 
II. 6. 6.  This  is  done  by  multiplying  II. 5. 5  by  a  factor  which  is  the 
ratio  of  the  mode  power  propagating  between  y  =  0  and  y  =  -t  (the 
region  where  As^(xy)  ^  0)  and  the  total  mode  power.  This  factor  can 
be  roughly  approximated  by 


(t  +  sin  hyt/py) 


(II. 6. 8) 


where  —  »  i-  the  penetration  depths  of  the  modes  into  the 

py  S 

upper  and  lower  boundaries  satisfy.  The  expression  for  the  coupling 


coefficient  between  two  channel  guides  thus  becomes: 


-kk- 


*10  =  %  -  -*  -  -1 


?  -P  8 

2a  p  t  e  * 
X  x 


(t  +  sin  hyt/py)(v  +  £-)(h*+v*) 


(11.6.9) 


For  a  very  veil  confined  mode 


t  »  —  and  w  »  — 
P  P 

y  x 


and  11.6.9  reduces  to 


*10  -  *01  -  -iK 


«  -1 


-p  S 

2h2pe  * 
x*x 

@w(h^p2) 


(II. 6. 10) 


which,  aside  from  its  negative  sign,  is  the  expression  derived  by 
Marcatili'2^ . 

Before  comparing  the  theoretical  expressions  for  the  coupling 
coefficient  with  the  experimental  results,  let  us  consider  a  geometry 
which  allows  an  accurate  measurement  of  the  coupling  coefficient. 


II. 7  Multichannel  Directional  Coupler  -  Coupling  Coefficient  Measurement 
In  Section  II. 3  we  have  considered  a  directional  coupler  made 
of  two  channels  in  which  the  power  oscillates  back  and  forth  between 
the  channels  as  a  function  of  propagation  distance.  We  consider  now  an 
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infinite  number  of  identical  guides  each  coupled  to  its  adjacent 
guides,  as  shown  on  the  left  of  Fig.  II. 1.  When  light  is  coupled  into 
a  central  channel  (defined  as  n  «*  o)  it  will  leak  to  the  sides. 

Because  of  the  large  number  of  channels,  it  is  clear  that  the  full  power 
will  not  appear  again  in  the  central  channel,  but  will  constantly  flow 
to  the  sides. 

The  coupled  mode  equation  for  the  nth  cnannel  is  given  by: 

dA  (z) 

nil  =  -i0An(z)  =  ikAn-1(z)  -  ikA^z)  (II. 7 .1) 

uZ 

where  B  is  the  propagation  constant  which  includes  tho  guide  attenuation  a 

P«Pr-if  (11.7.2) 

n  is  the  guide  number  (n=0,  ±1,  ±2, . . . ),  and  K  (a  real  number)  is 
the  coupling  coefficient  between  two  adjacent  guides.  (The  coupling 
coefficient  between  non-ad jacent  guides  is  negligibly  small).  When 
light  is  coupled  into  the  0  channel  only,  the  boundary  conditions  become 

A0(0)  -  1  V>(0)  • 0  <rI-T-3) 

The  solution  of  II. 7.1  is 

An(z)  =  (-i)n  Jn(2Kz)  e‘iPz 


(11.7,4) 


—1*6— 

and  the  power  flow  in  the  guides  therefore  is: 

PB(*>  =  Vz)  An(z)  =  Jn(2Kz)  e*°B  (11.7.5) 

where  Jn  represents  the  Bessel  function  of  nth  order. 

A  multichannel  direction  coupler  was  fabricated  by  proton 
implantation  through  an  appropriate  gold  mask  on  the  surface  of  the 
GaAs  substrate  as  shown  in  Fig.  II. 4b.  The  width  of  the  embedded  channels 
can  be  estimated  by  the  width  of  the  notches  in  the  mask  to  be  about 
2.^1  separated  by  3*9n  •  Their  depth  is  determined  by  the  known 
penetration  of  the  300  KeV  protons  in  GaAs  which  is  3p  .  The  index 
discontinuity  as  measured  in  a  fashion  described  in  Chapter  IV,  is 
about  .005  at  1.15u  and  about  10)6  larger  at  1.06u  . 

Figure  II. 1  describes  the  multichannel  directional  coupler 
with  the  input  light  coupled  into  a  central  guide  and  the  output  light 
scanned  at  different  sample  lengths.  The  argument  of  the  Bessel 
functions  at,  each  of  the  scans  was  determined  by  fitting  the  square 
root  of  the  normalized  intensities  of  the  light  in  the  different  channels, 
to  a  plot  of  the  Bessel  functions  as  shown  in  Fig.  II. 10.  The  experiment 
was  performed  with  light  having  a  wavelength  of  \  =  1.15u  (from  a 
HeNe  laser)  and  \  -  1.06p  (a  Nd:YAG  laser).  The  coupling  coefficient 
can  be  deduced  quite  readily  and  accurately  now,  simply  by  plotting  the 
argument  of  the  Bessel  functions  against  the  distance  of  propagation, 
and  measuring  the  slope. 


m  - 
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Fig.  II. 10  Sketch  of  various  orders  Bessel  functions  to  which  the 
square  root  of  the  amplitudes  in  Fig.  II  1  are  fitted. 
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Using  11.7.4  the  coupling  coefficient  is  thus  given  by: 

K  =  slope/2  (11.7.6) 

A  series  of  measurements  at  1.15u  and  1.06u  is  plotted  in  Fig.  II. 11. 
By  measuring  the  slopes  we  find  that  for  the  same  directional  coupler, 

K1.15u  =  °*52  ±  °'01  ““  1  h. 06u  =  °*3°  *  °*01  n,m"1  (i;e-7-T) 

It  has  also  been  found  that  different  polarizations  (E*  or  E^)  of  the 
input  beam  had  no  noticeable  effect  on  K  . 

Only  a  qualitative  comparison  between  theory  and  experiment 
can  be  made  with  these  samples.  The  modes  in  these  guides  are  not  well 
confined  modes,  to  which  the  theory  of  Section  II. 6  pertains.  Besides, 
the  dielectric  discontinuity  may  not  be  uniform  across  the  guide's 
cross  section  and  the  accuracy  in  measuring  it  is  not  high.  However, 
the  coupling  coefficients  calculated  on  the  basis  of  the  dimensions 
given  earlier  in  (II. 6. 9)  are  quite  close  to  the  actual  values.  The 
theoretical  results  are  (for  both  polarizations): 

K1.15  =  °-32  ^.06  =  °*21  nm’1  (11.7.8) 

Considering  the  limitations  mentioned  above,  these  results  are 
in  good  agreement  with  (II .7 .7). 

The  absorption  of  the  propagating  light  is  relatively  high  in 
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these  ion  implanted  waveguides.  It  amounts  to  about  5  cm”  which  is 
the  same  as  the  absorption  of  planar  ion  implanted  guides'  ' .  It  is 
caused  by  the  absorption  of  the  defect  centers  created  by  the.  implanta¬ 
tion  and  by  the  free  carrier  absorption  in  the  substrate.  Other 
methods  for  fabrication  of  directional  couplers  will  be  discussed  in 
Section  II. 9. 

II. 8  The  Sign  of  the  Coupling  Coefficient 

Let  us  consider  again  the  dual  directional  coupler  shown  in 
Fig.  II. 3.  The  amplitude  of  the  fields  in  the  two  guides  are  given  by 
II. 3. 5  and  11.3.6  as: 

AQ(z)  =  cos  Kz  eipz  (II. 3-5) 

Ax(z)  =  -i  sin  Kz  e^z  (11.3.6) 

From  these  equations  we  learn  that  a  distinct  phase  difference  exists 

between  the  two  guides.  We  also  find  that  the  phase  of  the  guide  in 

which  the  power  is  increasing  will  always  lag  90°  behind  the  phase  of 

the  guide  in  which  the  power  is  decreasing.  Formally  this  is  a  direct 

correspondence  of  the  sign  of  the  coupling  as  given  in  11.4.12. 

Hiysically  the  reason  for  the  time  lag  is  the  necessary  phase  relation 

between  the  polarization  (caused  by  the  field  in  guide  zero),  and  the 

field  in  guide  1,  if  power  is  to  be  generated  in  guide  1.  It  is  well 
(9) 

known' that  power  dissipation  in  a  dielectric  occurs  when  the  polari¬ 
zation  lags  the  field.  Therefore,  in  our  case,  feu*  power  generation 
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in  guide  1  the  field  there  has  to  lag  the  polarization,  which  is  caused 
by  the  field  of  guide  0  and  has  its  phase.  This  reasoning  thus  gives 
the  same  result  as  the  above  equations. 

The  phase  relation  between  the  two  guides  might  be  of 
importance  in  an  optical  circuit  in  which  the  phase  of  the  light  as 
well  as  its  amplitude  may  be  used  to  carry  information.  In  Ref.  2 
the  coupling  coefficient  appears  to  have  an  opposite  sign  to  the  one 
in  11.4.12.  Such  a  sign  will  cause  a  time  lead  in  guide  1  rather  than 
time  lag.  The  following  experiment  was  carried  out  to  find  whether  the 
phase  difference  is  time  lead  or  time  lag. 

Phase  information  can  be  acquired  by  interference.  In  this 
case  the  interference  was  performed  in  the  fashion  described  in  Fig.  11.12. 
Drawing  (b)  describes  a  geometry  in  which  the  input  beam  is  focused 
with  an  objective  into  the  central  channel  of  the  multichannel  directional 
coupler.  The  input  face  at  the  channels  is  in  the  focal  plane  of  the 
objective  and  the  spot  size  of  the  beam  is  small  enough  to  couple  into 
one  channel  only.  The  distribution  of  light  at  the  out-out  plane, 
caused  by  the  coupling  effect,  is  shown  at  the  top  of  the  drawing. 

By  moving  the  objective  in  or  out,  the  focal  plane  will  not 
coincide  with  the  input  plane  of  the  sample,  the  spot  size  will  increase 
and  light  will  be  coupled  into  the  adjacent  guides  as  well  as  the  central 
guide.  When  the  lens  is  moved  toward  the  sample,  the  curvature  of  the 
phase  fronts  of  the  beam  is  such  that  the  light  coupled  into  the  side 
channels  at  the  input  plane  has  a  time  lead  (as  shown  in  drawing  a). 
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On  the  other  hand,  when  the  beam  is  moved  away  from  the  sample  (drawing 
c)  the  light  coupled  into  the  adjacent  channels  suffers  a  time  lag. 

This  input  light  '•111  interfere,  as  it  propagates  with  the  light 
coupled  into  the  channels  from  the  centred  one.  Ike  intensity  of  the 
light  in  the  output  plane,  with  respect  to  the  central  guide  indicates 
where  this  interference  was  a  constructive  or  destructive  one.  For 
example,  in  Fig.  II.  12a  we  find  that  the  intensity  of  the  light  in  the 
adjacent  channels  is  close  to  zero,  while  in  (c)  it  is  very  strong. 

This  indicates  that  there  is  a  constructive  interference  in  case  (c) 
and  the  light  coupled  via  the  directional  coupling  effect  suffers  a 
time  lag. 

II. 9  Ridged  Channel  Waveguides  and  Directional  Couplers 

A  variety  of  thin  films  suitable  for  optical  waveguiding  has 

been  reported  so  far.  These  films  include  epitaxial  layers  of  high 

resistivity  on  low  resistivity  GaAs^10^,  composite  structure  of  GaAs- 

GaAlAs^11^  composite  structure  ADP-KDP^12^,  and  single  crystal  garnet 
(13) 

films'  .  !Ehese  layers  can  be  grown  with  a  high  degree  of  purity  and 
offer  the  attractive  feature  of  modulation^^,^',^t|^ . 

All  these  layers  are  planar  guides  which  confine  the  light  in 
one  dimension  only,  and  therefore  as  such  are  not  compatible  with  the 
concept  of  optical  circuitry.  To  form  a  channel  waveguide  in  which  the 
radiation  is  confined  in  the  two  dimensions  perpendicular  to  the 
direction  of  propagation,  it  is  possible  as  shown  earlier,  to  use  ion 
implanatation  or  diffusion .  It  is  interesting  to  investigate  the 
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possibility  of  fabricating  a  channel  guide  by  a  removal  of  the  super¬ 
fluous  sections  of  the  epilnyer  (as  shown  in  Fig.  11.13).  Because  of 
the  large  dielectric  discontinuity  the  scattering  losses  in  this  guide 
are  very  sensitive  to  the  smoothness  of  its  walls  as  well  as  its  top 
surface.  For  example,  a  3u  wide  waveguide  in  GaAs  will  suffer  roughly 
a  loss  of  4  cm”1  when  the  roughness  of  the  walls  is  50CA  rms  . 

This  rough  estimate  based  on  the  Rayleigh  criterion  indicates  that 
special  care  has  to  be  taken  in  the  fabrication  of  these  guides.  An 
extreme  case  of  rough  walls  and  top  surface  is  shown  in  Fig.  11.14a. 

This  guide  (3u  high,  7(i  wide)  scattered  away  most  of  the  propagating 
light  in  a  sample  less  than  1  am  long. 

Smoother  channels  were  fabricated  by  exposing  the  photoresist 
with  a  better  mask,  and  ion  polishing  the  sample  after  the  ion  machining 
(details  are  given  in  Chapter  IV).  Figure  II.l4b  shows  a  channel 
( 1.4n  high,  2  (i  wide)  with  much  smoother  walls. 

The  number  of  modes  supported  by  a  ridged  guide,  as  discussed 
in  Section  II. 6,  can  be  controlled  for  a  given  dimension  by  choosing 
an  appropriate  guide  substrate  dielectric  discontinuity  or  by  covering 
the  guide  with  an  appropriate  index  material. 

We  have  discussed  earlier  directional  couplers  made  of  single 
mode  waveguides  embedded  in  the  surface  of  GaAs.  In  that  case  the 
coupling  between  two  adjacent  guides  is  caused  by  the  overlap  of  the 
propagating  modes  as  shown  in  Fig.  II. 15a.  Figure  II. 15b  shows  single 
mode  ridged  waveguides  in  which  the  modes  are  very  well  confined  in  the 
x  direction  (becarse  of  the  large  dielectric  discontinuity)  and  two 
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#13  Channel  waveguide  fabrication  by  removal  of 
superfluous  sections  of  an  epilayer. 


(a)  An  extreme  example  of  a  guide  with  rough  walls  and 
top  surface.  (The  guide  is  3(i  high  and  7p  wide.) 

(b)  Good  quality  channel  waveguide  fabricated  by  using 
holographically  prepared  mask  and  ion  polishing. 

(The  guide  is  1.4p  high  and  2 p  wide.) 


y 


(a)  (b)  (c) 


Fig.  11.15  Directional  couplers  fabricated  in  GaAs. 

(a)  Single  mode  embedded  guides  directional  coupler. 
>1  Two  single  mode  ridged  guides  with  no  coupling, 
(c)  Ridges,  guides  with  increased  degree  of  coupling. 
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close  ly  spaced  guides  will  thus  have  a  negligible  coupling.  To  increase 
the  coupling  between  the  guides,  only  a  partial  removal  of  the  epilayer 
between  them  was  performed  as  shown  in  Fig.  II.  15c.  Fig.  11.16  is  an 
indication  of  the  coupling.  It  shows  the  cross  section  of  a  high 
resistivity  GaAs  epilayer  (about  8n  thick)  that  was  machined  down  to 
the  substrate  to  £>rm  a  large  number  of  isolated  channels.  Light 
coupled  into  one  channel  emerges  at  the  other  end  of  the  sample  from 
that  channel  only  (bottom  photograph).  However,  when  the  same  epilayer 
is  only  partially  machined  (Fig.  II.l6b),  light  coupled  into  one  guide 
emerges  from  three  guides. 

II. 10  Directional  Coupler-- Switch  Modulator 

As  mentioned  earlier,  one  of  the  most  exciting  applications  of 
the  directional  coupler  is  the  switch  modulator  in  which  the  amount  of 
light  coupled  from  one  channel  to  the  other  is  controlled  by  electric 
signal  applied  to  the  guides.  Such  a  device  can  be  used  to  modulate 
the  light  propagating  in  one  guide  or  to  switch  light  from  one  channel 
to  another.  Ihe  latter  property  is  an  important  mean  for  multiplexing 
(or  demultiplexing)  two  signals  into  one,  or  for  switching  light  in 
an  opto-electronic  switchboard  of  an  optical  communication  network. 

The  merit  of  a  modulator  is  determined  by  two  factors:  the 
sptsd  at  which  it  can  impose  the  information  on  the  light  beam  and 
the  electrical  power  it  consumes  in  order  to  do  so.  The  speed  of  the 
crystalline  electrooptic  modulator  is  determined  by  the  response  of  the 
electrical  circuit,  and  this  response  is  limited  by  the  capacitance  of 
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MULTIMODE  DIRECTIONAL  COUPIER 


Pig.  11.16  Directional  coupler  fabricated  by  ion  machining  in 


GaAs  epilayer. 

(a)  Isolated  channels. 


(b)  Coupled  channels 
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the  device.  The  capacitance  of  a  3  x  3n  rectangular  guide,  1  cm 
long,  is  about  «  10  pF  allowing  sub nano- second  response.  We  shall 
thus  limit  our  attention  to  the  relative  magnitude  of  voltage  necessary 
for  deriving  different  schemes  of  modulation.  First  we  recall  that  for 
the  linear  electrooptic  effect 

fin  «  E  =»  6n  «  V  (II. 10.1) 

where  fin  is  the  index  of  refraction  caused  by  the  applied  field  E 

and  where  V  is  the  applied  voltage.  We  also  recall  that  for 

(9) 

polarization  modulation v  (in  which  the  polarization  of  the  propagating 

O 

light  is  rotated  by  90  as  a  result  of  the  applied  field)  the  change 
60  in  the  propagation  constants  of  the  TE  and  TM  modes  must  be  such 
that: 

(6PTE  '  6fWL  =  *  (II. 10.2) 

where  L  is  the  length  of  the  modulator. 

For  the  sake  of  simplicity  we  assume 

=  0  (II. 11.3) 

a 

and  that  fip^  can  be  written  as 

60^  =  6n.k0  (II.  10.4) 

Using  II. 10. 4  and  II. 10. 2  we  can  write  that: 


iridi 
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6n  =  —T 


ko^ 


(II. 10. 5) 


Our  task  mill,  compare  this  6n  which  is  required  for  a  loo* 
modulation  In  the  pollination  modulator  (an  well  as  mode  conversion 
modulator'  ?>)  with  the  one  needed  for  a  direction^  coupler  modulator. 

Imt  u.  therefore  consider  a  typical  example.  It  is  the 
directional  coupler  whose  cross  section  is  shown  in  Fig.  n.17.  u* 
material  is  assumed  to  he  <WAs  (n  .  3.5).  Because  of  the  relatively 
large  dielectric  discontinuity  (in  =  .01)  the  modes  are  well  confined, 
Fielding  (with  the  use  of  11.6.9  for  the  3  x  3U  guides)  the  following’ 
value  for  the  coupling  coefficient  at  X  =  i.i5u 


K  =  0.155nm~^ 


(II. 10. 6) 


We  choose  the  modulator  length  to  he  the  length  necessary  for  complete 
transfer  of  power.  According  to  II. 3. 9 


^  z  lOmn 


(II. 10.7) 


urns  light  entering  the  device  at  the  left  channel  will  emerge  from  the 
one  on  the  right.  In  order  to  turn  the  device  into  a  switch  modulator 
Marcallti  2  has  suggested  changing  the  coupling  coefficient  (by  changing 
the  dielectric  disccntinulty,  through  the  electrooptic  effect)  to  a 
new  value  such  that 


k  ~  n 
K1  E  L 


(II. lo. 8) 
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Accarding  to  II. 3. 7,  this  will  cause  the  light  to  be  coupled  back  to 
the  original  channel  resulting  in  a  10O£  switch  modulation.  The  coupling 
coefficient  can  be  varied  by  applying  field  to  the  section  between  the 
two  guides  (schematically  via  electrode  b  in  Fig.  II. 17),  or  alternatively 
by  applying  the  field  to  the  guides  (electrodes  a  and  c).  To  obtain 
the  above  value  of  the  dielectric  discontinuity  of  the  guide  has 
to  be  reduced  from  An  =  .01  to  An  =  .007.  Thus  the  change  in  the 
index  of  refraction  required  for  switching  is: 


6n  b  .003  =  52  ~r  (II.10.9) 

K0  h 

The  conclusion  drawn  from  this  extreme  example  is  very  disap¬ 
pointing.  It  indicates  that  the  electrical  power  necessary  to  drive 
this  directional  coupler  modulator  is  2500  times  larger  than  the  driving 
power  of  a  regular  polarization  modulator  II. 10. 5.  A  smaller  6n  can 
be  achieved  far  a  given  length  geometry  by  decreasing  the  dielectric 
discontinuity  and  increasing  the  separation  between  the  guides.  This 
however  reduces  the  confinement  of  the  propagating  modes. 

In  order  to  devise  a  new  scheme  that  will  allow  directional 
coupling  switching  with  less  drive  power,  let  us  consider  again  the 
coupled  mode  equations  (II. 3. 2  and  II. 3. 3),  but  this  time  allowing  for 
different  propagation  constants  and  3^  for  the  two  guides. 

dA  (z) 

0  =  -iP^s)  -  iK  Ax(z)  (II. 10. 10) 


PteceM  H* wa* 
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dA.(2) 

— -  =  -10^(2)  -  iKAQ(2)  (II. 10. 11) 


Using  the  regular  boundary  conditions 


A0(0)  =  1  Ax(0)  =  0  (11.10.12) 


This  solution  becomes  identical  to  11.3*7  and  11.3.8  when  Ap  =  0  . 

Going  back  to  our  previous  example,  where  with  two  identical 
guides  (AP  =  0)  we  had  a  complete  transfer  of  power  in  a  length  L 
(A^  (L)  =  1),  it  is  possible  by  applying  field  to  one  of  the  guides 
only  to  destroy  the  identity  between  the  guides  so  that  Ap  ^  0  any 
more.  According  to  II. 10.14  if 


(II. 10.16) 
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we  will  get 

AX(L)  =  0 
Aq(L)  =  1 

Ibis  means  that  for  such  a  A3  there  is  no  power  transfer,  and  thus  a 
10056  switch  modulation  is  achieved.  Equations  II.10.l6  and  11.10.7 
give  the  value  of  A3  needed  for  destroying  the  coupling  between  the 
guides ; 

A3L  =  «  (II.lO.i7) 


and  the  necessary  change  6n  is  thus 


6n  »  1.73  ITT  =  i-10"4 

V 


(II.10.3B) 


This  6n  is  much  smaller  than  the  previous  one  II.  10. 9  and  is  approaching 
the  value  required  for  polarization  modulation.  It  is  interesting  to 
note  that  II.10.l8  applies  to  the  ridge  directional  couplers  described 
in  Section  II. 9  as  well  as  the  embedded  ones. 

A  directional  coupler  can  be  made  polarization  sensitive  by 
choosing  an  appropriate  guiding  plane  (for  example  100  in  GaAs)  in 
which  application  of  field  causes  a  certain  6n  for  TE  modes  but 
6n  =  0  for  IM  mode  (or  vice  versa).  Thus  in  our  example,  before  the 
application  of  the  field,  TE  as  well  as  IM  modes  are  coupled  from  one 
guide  to  the  other.  However  with  the  application  of  the  field,  the 
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coupling  of  the  IE  mode  is  destroyed  while  the  coupling  of  the  3M  mode 
remains  the  same.  Such  a  device  with  a  D.C.  electric  field  is  capable 
of  separating  the  two  polarizations  of  an  incoming  optical  signal  into 
two  different  channels.  As  such  it  can  turn  a  polarization  modulator 
(or  a  mode  converter)  into  a  directional  coupler  switch,  as  shown  in 
Fig.  11.16a.  Ibis  composite  device  will  have  a  low  modulating  power 
together  with  a  directional  coupler  switch  capability.  Reversing  the 
pulse’s  direction  of  propagation  (Fig.  II.lBb)  enables  the  devices  to 
combine  (multiplex)  into  one  output  channel  pulse,  with  the  appropriate 
polarizations  from  two  input  channels. 

We  can  apply  electric  fields  in  a  semiconductor^1^  guide  by 

back- biasing  a  metal-semiconductor  junction  (Schottky  barrier).  The 

field  is  supported  in  a  depletion  region  which  is  swept  free  of  carriers 

to  a  depth  d  ,  given  by  d  =  cE/eN,  where  N  is  the  free  carrier 

concentration  (assumed  independent  of  depth)  and  E  is  the  peak  value 

of  the  electric  field.  In  this  depletion  region  the  field  falls  off 

linearly  with  depth.  The  maximum  depth  to  which  a  field  can  be  applied 

in  such  a  junction  is  given  by  the  field  at  which  the  material  breaks 

down.  In  GaAs  the  breakdown  field  is  roughly  independent  of  doping 

5 

level,  and  has  a  value  Eg  5  x  10  V/cm.  This  means  that  the  maximum 
depletion  width  is  inversely  proportional  to  the  carrier  concentration. 
Numerically  the  maximum  depth  dffi  in  microns  is  related  to  the  carrier 

-3 

concentration  N  in  cm  by 

dm  =  3.5  x  1016/N  (II. 10.19) 
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Fig.  11.18  (a)  Combination  of  a  phase  modulator  and  a  polarization 

sensitive  directional  coupler. 

(b)  Polarization  sensitive  directional  coupler  as  a 
multiplexer. 
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This  indicates  that  in  order  to  penetrate  about  3u  into  the  guiding 
layer,  its  concentration  must  he  lower  than 

N  <  10l6cm" 3 

Unfortunately  the  free  carrier  concentration  of  the  proton  implanted 
waveguides  after  the  anneal  (which  is  necessary  for  low  optical 
attenuation)  is  about  10  'em  ,  thus  preventing  the  penetration  of 
the  field  into  the  bulk  of  the  guide.  GaAs  epitaxial  layers  with  low 
carrier  concentration  have  been  grown ^ 10  \  but  are  not  yet  readily 
available.  A  ridged  directional  coupler  (Section  II. 9)  made  of  such  an 
epitaxial  layer  may  be  a  suitable  candidate  for  a  modulator. 

In  a  GaAs  guide  in  the  (..00)  plane  1M  modes  are  unchanged  by 
the  applied  field.  TE  modes,  however,  experience  a  change  in  refractive 
index 

6n  =  &  n3  r4lE  ,  (II. 10. 20) 

where  r4l  is  the  electrooptic  coefficient  and  n  is  the  index  of 
refraction.  Substituting  the  value  of  n^  and  r  we  get: 

6n  =  3*10"9-E  (11.10.21) 

where  E  is  the  applied  field  in  V/cm. 

-4 

To  achieve  a  6n  =  10  which  is  necessary  for  the  directional 
coupler  switching,  we  need  an  average  electric  field  of: 


-70- 


li  v  v 

E  =  3.3*10  ^  =  3.3  zr 

cm  |i 

and  voltage  required  for  a  3u  thick  guide  is  thus 

V  =  10V 

In  the  next  section  we  shall  discuss  the  use  of  this  switch 
and  a  similar  device  for  multiplexing  and  demultiplexing  of  optical 
pulses. 

II. 11  Light  Multiplexing  by  Directional  Coupling 

The  high  optical  frequency  of  the  laser  radiation  enables  it 
to  carry  large  amounts  of  information.  However,  imposing  a  high  rate 
of  information  on  a  laser  beam  is  a  difficult  matter.  One  of  the  ways 
around  this  difficulty  is  to  modulate  the  information  on  a  number  of 
beams  at  a  relatively  moderate  rate,  and  then  to  multiplex  these  beams 
into  one  beam  as  shown  in  Fig.  II.  19a,  We  shall  assume  for  the  following 
discussion  that  the  laser  radiation  is  in  the  form  of  very  short  pulses 
with  a  period  T  ,  and  that  they  appear  in  the  different  channels  with 
the  appropriate  shift  in  time  so  that  they  can  be  multiplexed  with  no 
overlap  between  them.  Such  pulses  can  be  generated,  for  example,  by 
mode  locking'1'7' an  Nd:YAG  laser,  and  the  appropriate  shift  in  time  between 
the  channels  can  be  achieved  by  coupling  the  pulses  to  the  channel  via 
optical  fibers  with  different  lengths.  Another  method  is  to  have  an 
array  of  injection  lasers  on  the  same  substrate  as  part  of  an  integrated 
optical  circuit,  where  each  one  of  these  lasers  is  pulsed  at  a  moderate 
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rate  but  with  an  appropriate  shift  in  time. 

Let  us  examine  the  operation  of  the  directional 
coupler  switch  described  in  the  previous  section  as  a  two  channel 
multiplexer.  Figure  II. 19b  shows  the  channels  into  which  the  two 
already  modulated  (pulse  code  modulation- PCM)  beams  are  coupled.  In 
order  to  combine  the  two  beams  into  one  that  will  propagate  in  the  lower 
guide,  we  have  to  see  to  it  that  when  a  pulse  appears  in  the  upper 
guide  it  is  completely  coupled  into  the  lower  guide,  but  a  T/2  seconds 
later  when  a  pulse  appears  in  the  lower  guide  it  is  not  coupled  into 
the  upper  one  but  keeps  propagating  in  the  lower  guide.  This  can  be 
accomplished  by  choosing  the  coupling  coefficient  and  the  length  of 
the  two  guides  to  be  such  that  with  no  application  of  electric  field  a 
complete  transfer  of  power  from  one  channel  to  the  other  is  possible 
(KL  =  |)  .  On  the  other  hand  to  prevent  the  coupling  of  pulses  from 
the  lower  to  the  upper  channel,  we  have  to  apply  field  with  the  appro¬ 
priate  strength  for  destroying  the  coupling  whenever  a  pulse  in  the  lower 
guide  is  due.  This  can  be  done  by  driving  the  directional  coupler  with 
a  sinusoidal  wave  whose  period  is  exactly  T  as  shown  in  Fig.  II. 19b. 

To  perform  multiplexing  the  same  configuration  can  be  used  but  the 
direction  of  propagation  of  the  pulses  is  reversed. 


Fig.  11.19  (a)  Multiplexing  scheme  for  a  higher  rate  of 


information  modulation. 


(b)  Directional  coupler  multiplexer. 
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III.  PERIODIC  CORRUGATIONS  IN  OPTICAL  WAVEGUIDES 
III.l  Introduction 

The  interaction  of  radiation  with  man-made  periodic  structures 
ia  of  great  importance  in  the  microwave  and  optical  regions.  The 
traveling  wave  tube  is  a  good  example  of  a  device  in  the  microwave 
region  which  is  based  on  this  kind  of  interactions.  In  the  optical 
region  one  finds  devices  such  as  the  diffraction  grating,  the  acousto- 
optic  modulator  and  deflector,  and  the  hologram. 

In  this  chapter  we  shall  discuss  the  interaction  between  periodic 
corrugations  on  the  surface  of  a  dielectric  waveguide  and  the  wave 
propagating  in  the  guide.  In  Chapter  V  we  shall  describe  methods 
which  were  devised  to  enable  us  to  fabricate  these  periodic  structures, 
with  periods  as  low  as  O.II4  (HOOA). 

In  general  a  wave  propagating  in  a  periodic  structure  consists 
of  space  harmonics.  These  space  harmonics  have  the  same  frequency  as 
the  original  wave  but  a  different  propagation  constant  p  .  Considering 
the  first  order  space  harmonics  only,  we  can  write 

Pspace  harmonic  =  P  ±  T  (III. 1.1) 

where  A  is  the  corrugation  period.  The  space  harmonics  are  thus 
shifted  in  k  space  by  the  length  of  the  k  vector  of  the  structure  -T1 


(111.1,2) 


If  the  propagation  constant  of  the  space  harmonic  happens  to  coincide 
with  the  propagation  constant  of  a  different  mode,  coupling  of  power 
from  one  mode  to  the  other  becomes  possible. 

Let  us  consider  the  guide  shown  in  Fig.  III. la;  the  u>-0 
diagram  is  shown  in  Fig.  III. lb.  The  guide  (which  is  discussed  in 
Appendix  I)  is  assumed  to  support  two  guided  modes  0  and  1.  The  propa¬ 
gation  constants  of  these  modes  at  a  given  cu  ,  are  0Q  and  0^  . 

As  can  be  seen  from  the  diagram  (we  have  assumed  n^  >  n^) 

Vo  >*0><il>  Vo  (in. 1.3) 

The  0's  are  less  than  rigk Q  in  order  to  allow  the  cosine  or  sine 
variation  of  the  modes  inside  the  guide  along  the  thickness  direction. 
Outside  the  waveguide  the  mode  profile  should  fall  exponentially  for  a 
guided  mode.  This  is  taken  care  of  by  requiring  0  >  nQk0  where  n^kQ 
is  the  free  wave  propagation  constant  in  the  substrate. 

In  this  waveguide  it  is  possible  to  cause  coupling  between  the 
two  guided  modes ^  by  having  a  corrugation  whose  period  satisfies 

=  7|  =  0Q  -  0X  (III. 1.4) 

This  condition  is  illustrated  on  the  right  branch  of  Fig.  III. lb  by 
T)^  .  The  value  of  the  coupling  coefficient  which  depends  on  the 
corrugation  depth,  the  modes'  profiles  and  so  on,  will  be  discussed 
later.  If  wo  decrease  the  corrugation  period  and  thus  increase  the 
value  of  its  vector  T)  ,  it  becomes  possible  to  couple  light  from  a 
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guide  d  mode,  say  the  0  mode,  to  &  continuum  of  radiation  modes  ^ . 

These  radiation  modes  have  a  propagation  constant  in  the  z  direction 
which  is 

Pr  <  ^  (III. 1.5) 

This  means  that  they  do  not  decay  exponentially  in  the  x  direction  and 
are  able  to  carry  power  away  from  the  guide.  The  case  is  illustrated 
by  T)2  in  Fig.  III. lb.  The  coupling  fircm  the  0  order  mode  will  be  to 
a  cluster  of  radiation  modes  whose  propagation  constant  in  the  z  direction 
is  close  enough  to 

Pr  =  0O  "  T)  (III. 1.6) 

This  effect  can  be  used  to  couple  light  in  and  out  of  a  waveguide  and 

(2 

is  known  as  the  grating  coupler'  *. 

If  A  is  further  decreased  so  that 

X  s  ^  =  2^0  (III. 1.7) 

coupling  occurs  between  a  forward  propagating  mode  and  the  same  backward 
mode.  This  effect  is  illustrated  by  T]^  in  Fig.  III. lb.  It  results 
in  a  distributed  mirror  whose  behavior  was  described  by  Fig.  1.5  and 
will  be  considered  later.  The  use  of  two  such  distributed  mirrors 
instead  of  two  cleaved  faces  may  help  in  extending  the  lifetime  of 
semiconductor  injection  lasers  in  the  case  where  catastrophic  mirror 
damage  is  the  dominant  failure  mode.  When  the  corrugation  extends 
along  the  whole  length  of  an  amplifying  medium,  a  distributed  feedback 
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laser  7  may  result.  In  these  lasers  there  are  no  discrete  mirrors 
and  this  makes  them  very  desirable  for  an  optical  integrated  circuit. 

In  this  chapter  we  shall  focus  our  attention  on  a  laser  geometry 
in  which  the  distributed  feedback  is  supplied  by  surface  corrugation, 
elaborate  on  its  modes,  and  conclude  with  a  description  of  the  first 
observation  of  such  a  laser  in  GaAs^. 


III. 2  Theory 

The  problem  of  coupling  by  surface  perturbation  between  HE 
modes  in  a  dielectric  waveguide  can  be  treated  in  the  same  fashion  as 
the  time  dependent  perturbation  in  quantum  mechanics.  The  solution  for 
the  perturbed  guide  can  be  expanded  in  terms  of  the  modes  of  the 
unperturbed  guide.  Assuming  that  the  zero  order  mode  enters  the 
perturbed  section  of  the  waveguide  at  z  =  0  (as  shown  in  Fig.  III. 2a), 
one  can  calculate  "transition  rates"  or  coupling  coefficients  to  dif¬ 
ferent  guided  modes,  to  radiation  modes  (modes  of  the  continuum)  and 
to  the  backward  going  mode.  This  approach  was  used  by  Marcuse ^  to 
calculate  mode  conversions  and  radiation  losses  in  waveguides. 

We  shall  apply  the  same  approach  used  in  Chapter  II  for  the 
calculation  of  the  coupling  coefficient  between  the  modes  of  two 
adjacent  guides.  We  start  by  writing  the  coupled  mode  equations  for 
forward  (+)  and  backward  (-)  going  zero  order  modes: 


dA^+)(z) 

dz 


(-%  a)  A^'kz)  -  iKe'iT]ZA^(z)  (III. 2.1) 
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Fig.  III. 2.  (a)  Distributed  mirror  of  length  L  . 

(b)  Definition  of  the  perturbation  in  the 

dielectric  constant  caused  by  the  corrugation. 


(III. 2. 2) 


»  (130  -  a)  A^(z)  +  iKe1"2  a£+' 
dz 

where  0Q  is  the  propagation  constant  of  the  mode,  a  is  the  gain  of 
the  medium,  and  K  is  the  coupling  coefficient  between  the  modes. 

K  appears  with  different  signs  in  the  two  equations  because  the  modes 
carry  power  in  opposite  directions  (see  Appendix  II).  In  addition, 

K  is  multiplied  by  a  phase  term  which  represents  the  spatial 

variation  of  the  coupling  perturbations.  As  we  shall  see  later,  T) 
is  exactly  the  k  vector  of  the  perturbation.  and  A^  are 

the  amplitudes  of  the  modes 

E^+)(y,z)  =  A^+)(z)e0(y)  (III.2.3) 

E^-)(y,z)  =  A^\z)£0(y)  (111.2.4) 

Hie  mode  profile  C^(y)  is  normalized  to  carry  one  unit  of  power,  so 
that  the  power  P^(z)  carried  by  E^(y,  z)  is  given  by 

p£+)(z)  =  |a£+)(z)|2  (III. 2. 5) 

The  rate  of  growth  of  P^(z)  because  of  the  coupling  can  be  derived 
from  III. 2.1  (we  assume  a  =  0  for  this  purpose) 


Rcfci'  \z)  [-i]fe‘iTV‘)(z)l‘}  (III. 2. 6) 


We  shall  derive  now  an  expression  for  the  power  growth  from  a  physical 
point  of  view,  and  by  a  comparison  with  111.2.6  find  expressions  for 
K  and  T\  . 


In  a  similar  fashion  to  what  we  did  in  Chapter  II,  we  use  the 
perturbations  in  the  polarization  of  the  backward  going  mode,  caused 
by  the  index  disturbance  at  the  surface,  to  drive  the  forward  going 
mode.  Coupling  can  take  place  only  between  TE  or  04  modes  but  not 
from  TE  to  1M,  because  the  driving  polarization  is  in  the  same  direction 
as  the  field  (we  assume  isotropic  media).  The  driving  polarization 
caused  by  a  TE  mode  is  given  by: 


[P^rt(y’Z't)]y  =  ^e(yfz)[E^(y,z,t)]y 


(III. 2.7) 


where  Ae(y,z)  is  the  perturbation  in  the  dielectric  constant  and  is 
given  in  Fig.  III. 2b.  Power  generation  for  the  (+)  mode  in  an 
infinitesimal  volume  dV  is  given  by: 


(III. 2.8) 


where  the  horizontal  bar  denotes  time  averaging.  Using  111.2.4  and 
time  dependence  of  e^  ,  111,2.8  becomes  after  tht  time  averaging: 


.(+) 


dV 


Re  A^(z)e0(y)l-injdE(y,z)eQ(y)A^(z)]* 


(III. 2.9) 


the  growth  of  the  total  power  P^'  ^  can  now  be  derived  by  simply 
integrating  over  x  and  y  .  We  obtain 


-83- 


dP<+)(z) 

dz 


OD 

=  $  Re  JAo+^z^"lu)|j  cJ(y)Ae(y,z)e0(y)dxdy|  Ao’^2)*| 


(III. 2. 10) 


Comparing  III.2.lOvith  111.2.6  we  find: 


Ke“i7Jz  =  £  jj  e2(y)Ae(y,z)dxdy  (III. 2. 11) 

•CD 

The  integral  in  111,2.11  can  be  simplified  if  the  perturbation  is 
enough  (a  «  X)  so  ve  can  approximate  eQ(y)  by  £Q(0)  .  We  ignore 
for  the  time  being  the  integral  over  x  because  we  consider  now  a 
planar  guide.  (The  x  integral  will  be  used  when  we  consider  a  channel 
laser.)  With  the  help  of  the  node’s  profile  expression  from  Appendix  I, 
we  get  after  simple  manipulations: 


where 


Ke 


-i7]z  _ 


h2  a 


(111.2.12) 


sgn[ argument]  = 


1  for  argument  >  0 
-1  for  argument  <0 


(III.2.13) 


Using  Fourier  expansion  we  can  rewrite  III. 2. 12  as 

Ke'iT|z  - £-5 — r[i  oos(^  z)  -  —  cos (3  ¥•)  +....1  (111.2.14) 

3”  J 

Only  one  term  from  the  left  side  of  111.2.14  will  give  synchronous 
contribution  when  substituted  into  the  coupled  modes  equations.  For 
first  order  coupling  from  to  A^  which  occurs  when  2p0 
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we  find 


«  2« 
T\  =T 


(m.2.15) 


and  the  coupling  coefficient  for  this  case  is 


K  = 


h2  a 

”f>0(t  +  5 +  1> 


(m.2.16) 


On  the  other  hand,  for  third  order  coupling  which  occurs  when 

2it 

3  -  20q  ,  the  synchronous  term  is  the  one  for  which 


ti  a  2  it 

<1  *  3  -j* 


(III. 2. 17) 


and  the  appropriate  coefficient  for  the  third  order  coupling  is 


K  = 


^u7T7|j 


(111.2.18) 


These  expressions  can  be  further  simplified  for  a  well  confined  zero 
order  mode  for  which 

h~!  (III.2.19) 


.  _  1  1 
t  »  -  +  - 
p  q 


(III. 2. 20) 


We  can  thus  rewrite  111.2.16  as 


«a 


K  »  s  t3 

P0  * 


(III. 2. 21) 
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An  additional  contribution  to  the  coupling  coefficient  may  arise  in  a 
case  where  the  gain  (or  the  loss)  a  of  the  medium  varies  periodically^. 
The  he  which  represent  gain  variations  is  imaginary,  and  the  combined 
Ae  becomes  complex: 

Ae(y,z)  =  AecorrU€ation(y,z)  +  i  ~  **(y,z)  (III.2.22) 

where  kQ  is  the  free  space  propagation  constant  and  Aa(y,z)  is  the 
variation  in  the  gain.  Using  this  expression  for  Ac(y,z)  one  has  to 
go  back  to  equation  III. 2. 11  and  evaluate  K  .  This  coupling  coefficient 
will  have  the  following  form 


K  =  K  ..  +  iK  . 

corrugation  gain 


(III. 2. 23) 


^corrugation  is  the  one  8iven  by  In-2-16- 
Now  that  we  have  obtained  the  value  of  K  and  1]  let  us 

return  to  the  coupled  mode  equations.  To  eliminate  from  the 

equations  we  redefine  Aq*^(z) 


(+)/  \  (0/  \ 

Aq  '(z)— ■  (z)e  0 


,<-)(,) - A(-)(a)eiCP°-iP) 


(111.2.24) 


where  A{J  is  the  phase  mismatch  constant 


A0  =  30  “  ^  T 


(III. 2. 25) 
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Substituting  all  that  into  III. 2. 2  and  III. 2. 2  we  get: 


-  =  (a-idP)A^  *(z)  -  iK  Aq  '(z) 


dz 


iA^(z) 

dz 


.(a-idpjA^z)  +  iKA<+)(z) 


(111.2.26) 


(III. 2. 27) 


This  pair  of  coupled  equations  will  be  our  starting  point  in  the  next 
sections.  We  recall  that  X  (for  a  well  confined  zero  order 
mode  and  a  «  \)  and  dp  are  given  by: 


K  * 


J 


ap  *  P0  x 


(111.2.28) 


III. 3  Distributed  Mirror 

We  consider  now  the  geometry  described  in  Fig.  III. 2a.  The 
(+)  zero  order  mode  propagating  in  the  guide  from  left  to  right 

2j f 

enters  at  z  =  0  a  corrugated  section  of  length  L.  If  2pQ  z 

reflection  occurs.  The  (+)  mode  decays  while  the  (-)  mode  grows 
as  indicated  in  Fig.  1.5.  Our  aim  now  is  to  find  how  strong  the 
reflection  is  as  a  function  of  the  guide  properties,  and  how 
selective  the  reflection  is  for  different  frequencies  that  propagate 
with  different  propagation  constant  -pQ  .  To  do  that  we  consider 
again  the  coupled  mode  equations  111.2.26  and  III. 2. 27,  taking  a  <=  0  . 
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dA^+)(») 

dz 

dA^"^(z) 

dz 


-1A0  A^+)(z)  -iKA<-)(z) 
iAp  A^(z)  +  iK  Aq+^(z) 


(III. 2. 26) 


(III. 2. 27) 


The  boundary  conditions  that  describe  the  case  under  consideration  are: 


Aq+)(0)  =  1 


A^_)(L)  =  0 


(III. 3.1) 


The  solution  is  given  by: 


A(-)(z)  =  -iK  sinh  y(L-z) 


Y  cosh  yL  +  iAp  sinh  yL 


(III. 3. 2) 


f.\  y  cosh  y(L-z)  +  iAp  sinh  y(L-z) 

Aq'(z)  - - 

y  cosh  yL  +  iAp  sinh  yL 


(III. 3. 3) 


where 


v  =  (K2  -  (4fS)2)* 


(III. 3.4) 


Under  phase  matching  conditions  (Zip  =  0)  the  Bragg  condition 


no  _  jiH 

2po  ‘  T 


(ill. 3. 5) 


is  fulfilled,  and  we  have  the  strongest  reflection: 


sinh  K(L-z) 

-i - 

cosh  KL 


cosh  K(L-z) 
cosh  KC> 


(111.3.6) 


(III. 3. 7) 
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A  plot  of  the  mode  powers  for  this  case  is  shown  in  Fig.  1.5.  The 
curves  of  |Aq"^|2  and  |A<->|2  are  plotted  for  KL  =  1.84  . 
According  to  III. 2. 21,  it  represents  the  reflection  of  a  well  confined 
mode  in  a  guide  with  t  =  3u  ng  =  3.6  \  =  0.8ji  a  =  50QH  and 

L  =  9ban  . 

The  selectivity  or  the  filtering  properties  of  the  mirror  can 
be  described  in  terms  of  dp  .  We  recall 

dp  =  p0(<o)  -  *  (m.2.23) 

it  can  be  approximated  by 

dp  »  (<u  — •  -  ^  -r-)  =  da)  — —  (III.3.8) 

c  A  c 

Thus  its  value  represents  deviation  from  the  exact  frequency  for  whu.ih 
the  Bragg  condition  III. 3. 5  is  fulfilled  (uvra__)  .  Figure  III. 3a  des 
cribes  the  reflection  properties  of  a  mirror  with  a  fixed  K  and  the 
phase  shift  associated  with  it.  The  different  curves  represent 
different  lengths.  As  can  be  seen,  the  longer  the  length  is,  the 
narrower  the  reflection  band  is.  For  a  finite  mirror  length  the  peak 
reflection  is  less  than  unity,  but  the  peak  is  relatively  wide.  Its 
width  APrefl  is  given  approximately  by 

ifwi.  -  +  4>4  (m.3.9) 

L 

An  infinitely  long  mirror  has  a  flat  total  reflection  band  with  a  width 
dp  of  2K  .  This  total  reflection  band  is  the  stop  band  of  this 
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Fig,  III.3.  (a)  Reflection  intensity  and  phase  shift  from  a  distributed 

mirror  as  a  function  of  A0  for  different  values  of  KL  . 
(b)  A  stop  band  in  the  mode  dispersion  caused  by  the  corrugation 


periodic  structure,  and  is  analogous  to  Gragg  scattering  of  Bloch 


electron  wave  in  a  crystal,  from  one  edge  of  the  BriUouln  zone  to  the 

(6) 

other  by  the  crystal  periodicity  .  The  dispersion  (<a  -  p)  diagram 
of  the  waveguide  thus  has  a  stop  band  and  is  modified  at  its  vicinity 
as  shown  in  Fig.  III. 3b.  The  midgap  frequency  is  the  Bragg  frequency 
for  which 

Pit 

VW  -  *  T  (III. 3. 10) 


The  height  of  the  energy  gap  is  the  frequency  region  over  which  the 
propagation  constant  becomes  complex,  and  thus  the  propagation  is 
forbidden.  (The  wave  is  not  absorbed  but  reflected.)  The  gap  region 
is  given  by 

APgap  =  2K  (III. 3. 11) 


and  if  we  take 
approximately: 


Bq((u)  ,  the  non-perturbed  propagation  constant*  a r 


P0(w)  * 


u>  _ 

—  n. 

c  2 


(III. 3. 12) 


we  get: 


(III. 3. 13) 


In  the  gap  the  real  propagation  constant  is  fixed  and  is  equal  to 

£  ^  .  Outside  the  stop  band  the  propagation  constant  is  real,  but  is 

modified  somewhat  to  accommodate  the  gap,  as  shown  in  Fig.  III. 3b. 

(l) 

Two  Bragg  reflectors  can  form  a  laser  cavity' ' ' .  Two  corrugated 
mirrors  can,  for  example,  be  used  instead  of  cleaved  faces  to  form  a 
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semicooductor  laser  cavity  as  shown  in  Pig.  III. 4*.  ftie  reflection  of 
the  mirrors  can  be  controlled  by  changing  their  length  or  the  depth 
of  the  corrugation  to  supply  the  appropriate  feedback  for  oscillations 
even  for  a  low  gain  medium.  We  assume  that  the  mirror's  properties 
do  not  change  much  in  the  presence  of  a  low  gain,  and  that  its  length 
is  much  smaller  than  the  laser's  length  (a  different  case  where  the 
mirrors  extend  all  over  the  laser  is  treated  in  the  next  section). 


L  «  L, 
m  v 


(111.3.14) 


This  means  that  the  spectral  reflectivity  of  the  mirrors  is  much  wider 
than  the  spacing  between  the  longitudinal  modes  of  the  laser.  We 
recall  that  the  width  of  the  reflectivity  in  &  space  is  given  by 


HI. 3. 9 


«Wror  *  2(k2  +  7?>4 

m 


(III. 3. 15) 


For  a  mirror  reflectivity  which  is  close  to  a  uTisy  ,»e  can  assume  that 


~  2 


(111.3.16) 


Ihu8  a  lower  limit  for  A(3  .  is 

mirror 


^mirror  >  2  ^ 


(III. 3. 17) 


On  the  other  hand  the  spacing  of  the  longitudinal  modes  is  given 
approximately  by:  (we  ignore  dispersion  because  of  the  gain  in  the 
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Fig.  III. 4.  (a)  A  laser,  formed  by  two  distributed  mirrors. 

(b)  Longitudinal  mode  spectrum  of  the  above  laser, 
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medium,  and  dispersion  in  the  phase  of  the  reflected  beam  from  the 
mirror  which  is  shown  in  Fig.  III. 3a  and  can  be  ignored  for  «  L^) 


longitudinal  modes  “ 


(111.3.18) 


The  longitudinal  modes  distribution  is  shown  in  Fig.  III.4b. 

The  location  of  the  modes  with  respect  to  the  center  of  the  mirror 
reflectivity,  which  is  important  for  longitudinal  mode  discrimination, 
depends  on  the  exact  distance  between  the  two  mirrors.  According  to 
111.3.6  the  reflected  wave  from  the  mirrors  suffers  a  phase  shift  of 
-90*  .  This  means  that  for  a  longitudinal  mode  to  coincide  with  the 
peak  of  the  reflection,  the  distance  between  the  mirrors  has  to  be 
such  that  it  compensates  the  -90*  phase  shift  for  the  Bragg  frequency. 
This  requires  that 


\  =  (2m*l)  =  (zmfl)  | 


(III.3.I9) 


Another  way  to  state  this  condition  is  that  the  relative  phase  of  the 
two  corrugations  has  to  be  180*  .  This  requirement  may  be  found 
difficult  to  fulfill  especially  when  the  propagation  constant  under 
the  corrugated  section  may  not  be  exactly  the  same  as  the  one  in  the 
bulk  of  the  guide.  In  addition  the  flatness  of  the  reflectivity  curve 
close  to  the  center,  reduces  the  ability  to  discriminate  between 
close  longitudinal  modes. 

When  the  corrugation  is  extended  over  the  total  length  of  the 

(4) 

laser,  a  distributed  feedback  laser'  results,  and  a  better  control 
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of  the  longitudinal  nodes  becomes  possible. 


Ill, 4  Distributed  Feedback  laser 

We  consider  now  the  geometry  shown  in  Fig.  III. 5a.  The 
corrugation  extends  from  z  =  to  z  =  £l  and  we  assume  that  the 
waveguide  has  a  gain  of  a  .  In  the  previous  section  ve  have  discussed 
the  reflection  properties  of  a  corrugated  structure.  In  this  section 
we  shall  discuss  its  lasing  properties  and  inquire  about  the  modes  of 
oscillation,  their  frequencies  and  threshold  gains. 

Our  starting  point  is  the  coupled  modes  equation  derived  in 
Section  III. 2. 


dA£+)(z) 


(a-i^)A<+)(z)-  iKA^fz) 


(111.2.26) 


dA^Cz) 


-(a-idP)A^"\z)  +  iKA^(z) 


(III. 2. 27) 


The  corrugation  is  responsible  for  the  coupling  mechanism  and  the 
coupling  coefficient  was  found  to  be  (for  a  well  confined  zero  order 
mode  and  corrugation  depth  a  «  \) 


(III. 2. 21) 


where  t  is  the  thickness  of  the  guide.  We  recall  that  AfJ  is  the 
deviation  of  the  unperturbed  propagation  constant  of  the  guided  mode 
from  n/A  . 
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Ap  =  PQ(o>)  -  i  &  (III. 2. 25) 

Prom  this  point  on  we  follow  the  original  discussion  of  the 
distributed  feedback  laser  by  Kogelnik  and  Shank ^ .  We  start  by 
imposing  the  appropriate  boundary  conditions.  Since  this  is  an 
oscillator  we  assume  that  there  are  no  incoming  waves.  Thus  the 
boundary  conditions  are: 

a£+)(-£l)  =  o  ,  a£")(?l)  =  0  (III. 4.1) 

A  solution  with  such  boundary  conditions  exists  only  when  the  following 
relation  is  fulfilled^ 


a  -  iAP  =  y  coth(yL)  (111,4.2) 

where 

Y  =  [(a  -  mf  +  K2]*  (111.4.3) 

The  complex  equation  111.4.2  which  can  also  be  written  as 


-V  +  (a- 1&P)  g2YL  ^ 
y  +  (a-iAp) 


(III. 4. 4) 


determines  for  a  given  K  and  L  the  eigenvalues  of  AS  (the  shift 
in  the  modes  frequency)  and  +he  corresponding  a  (the  threshold  gain) 
for  which  oscillations  can  take  place.  The  solution  for  AQ(z)  is 
given  by 

Aq+  ^(z)  =  sinh  y(z  +  £l)  (111.4.5) 
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Aq‘^(z)  =  sinh  y(z  -  fcL) 


(III. 4.6) 


and  is  shown  in  Fig.  III.5b.  As  can  be  seen  in  the  fi_gure,  each  wave 
starts  with  a  zero  amplitude,  at  a  different  side  of  te-he  .aser,  The 
wave  grow 8  as  it  propagates  because  of  the  energy  cent-  inuously  fed 
into  it  by  the  other  wave  and  because  of  the  gain  whies li  exists  in  the 
medium. 


In  the  experimental  cases  to  be  discussed  latter  oc  »  K  so 
we  confine  our  attention  now  to  this  high  gain  situati.  on.  Thus,  y  , 
which  is  given  by  111.4.3,  can  be  expanded  as 


Y  a  (a  -  idp)  + 


K2 

2(a  -  idp) 


(m. 4.7) 


and  Equation  III. 4. 4  can  be  written  approximately 


4(a  -  idp)‘ 


=  1 


(III. 4.8) 


The  phases  of  this  equation  give  us  the  phase  constants  sund 
thus  the  frequencies  of  the  modes  of  oscillation  (by  qLxiantizing  AP)  , 
while  the  absolute  value  of  this  equation  gives  the  tlvrreshold  condi¬ 
tions  for  the  different  inodes  (by  solving  for  a)  .  TJtx  distinct 
values  of  (Ap)m  for  which  oscillations  are  possible  mare  thus  given  by: 

a^ap)  ^(ap)  L 

2(ap)  L  -  arctan  -  -  *  -= - =r  =  (2»L  )it  (111.4.9) 

m  ^-(AP)*  (&(AP  )l 

where  m  is  an  integer  which  denotes  the  longitudinal,  mode  number. 


-«•** - - -  ^ 


This  equation  is  not  in  a  closed  form  because  has  notvMen 
determined  yet.  However,  if  we  assume  that  K  ,  AfJ  «  a  we  immediately 


<‘f»B  - 


mu  +  \  it 


(III.4.10) 


Using  111.3.8  we  find  that  for  such  a  laser,  the  spacing  of  the 
longitudinal  modes  is  approximately  the  same  as  the  spacing  of  a  regular 
two-mivror  laser  of  length  L  .  We  have 


Frequency  spacing  of  longitudinal  modes 


(III. 4. 11) 


We  also  find  that  we  do  not  have  a  mode  for  which  A£  =  C  .  This  means 
that  the  laser  will  not  oscil'^te  with  a  frequency  which  is  exactly 
the  Bragg  frequency.  A  physical  explanation  of  this  effect  will  be 
discussed  later. 

The  threshold  gain  for  the  m  longitudinal  mode  can  be  derived 
from  III. 4. 8  by  taking  the  absolute  value  of  the  equation.  It  yields 


°m+  ^m 


(111.4.12) 


The  vaLue  of  the  gain  can  be  found  either  by  using  the  approximate  value 

of  (A{i)m  (111.4.9)  or,  more  accurately,  by  solving  111.4.12  in 

conjunction  with  111.4.9  for  (A{3)  and  a  .  The  broken  line  in 

mm 

Fig.  III. 6  describes  a  specific  example  of  a  distributed  feedback  laser 

-1 

of  length  L  =  1  mm  and  K  =  2  cm  .  It  shows  the  required  threshold 
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gain  and  the  deviation  of  different  longitudinal  modes  from  the  Bragg 
frequency.  Since  A0  is  not  much  smaller  than  a  the  modes  do  not 
coincide  with  the  c/2nL  s pacings,  but  are  further  pushed  to  the  sides. 
The  selectivity  of  the  device  is  evident  from  higher  threshold  gain 
needed  for  the  higher  order  modes. 

As  we  have  observed,  the  laser  does  not  have  a  mode  for  which 
Af5  =  0  .  This  can  be  explained  in  the  following  manner.  We  recall 
that  according  to  111.3.6  the  reflected  wave  from  a  distributed  mirror 
suffers  a  phase  shift  of  -90°  ,  when  its  frequency  is  exactly 
“teragg^  =  ®1US»  we  88  811  example  a  wave  propagating 

to  the  right,  the  fraction  of  it  which  is  reflected  backward  and  then 
reflected  again  f  rward  will  suffer  a  total  phase  shift  of  -180*  .  This 
means  that  it  will  interfere  destructively  with  the  original  wave,  and 
prevent  oscillation  at  that  particular  frequency. 

It  is  interesting  to  note  that  this  problem  does  not  exist  when 
the  coupling  mechanism  is  provided  by  periodic  gain  variation.  According 
to  III. 2. 22,  in  this  case  K  is  imaginary,  giving  rise  to  a  zero  phase 
shift,  and  oscillation  exactly  at  the  Bragg  frequency  is  allowed^. 

This  mode,  for  which  A(3  =  0 ,  would  clearly  have  the  lowest  threshold 
gain. 

We  have  noted  earlier  that  a  corrugated  waveguide  has  a  stopband 
of  frequencies  (Fig.  III. 3b),  in  which  propagation  is  not  allowed.  The 
gap  increases  with  increasing  K  .  When  it  becomes  comparable  to  the 
c/2nL  s pacings,  it  starts  pushing  the  longitudinal  modes  to  the  sides, 
preventing  oscillations  inside  the  stop  band.  A  small  amount  of  pushing 
is  evident  in  the  second  example  in  Fig,  III. 6.  The  solid  lines 


describe  the  modes  of  a  laser  with  K  =  6  cm"  *  and  L  =  1  mm  .  Also 
evident  are  the  lower  threshold  gains  and  the  better  selectivity  of 
the  laser. 

III. 5  Distributed  Feedback  Laser  in  GaAs  -  Experiment 

Fundamental  Bragg  coupling  of  a  waveguide  laser  requires  a 
corrugation  with  a  period  of  ~  where  XQ  is  the  free  space 

oscillation  wavelength  and  the  index  of  refraction  of  the  waveguide. 
For  a,  ~  3.6  and  X  =  .82|i  the  requisite  period  is  ~  O.II41  . 

A  technique  for  the  fabrication  of  such  small  period  gratings 
in  solid  substrates  was  developed.  It  consists  of  ion  milling  a 
photoresist  grating  into  the  substrate,  which  is  GaAs  in  this  case. 

The  grating  in  the  photoresist  is  formed  by  exposing  it  with  an  inter¬ 
ference  pattern  of  two  laser  beams.  In  order  to  achieve  the  small 
period  grating  U.V.  laser  (HeCd  Xq  =  3250A)  was  used  and  the  wavelength 
was  further  reduced  by  sending  the  beams  through  a  prism  (n  *“  1.5) 

(the  prism  is  attached  to  the  photoresist  with  index  machining  oil 
between  them).  A  scanning  electron  microscope  photograph  of  a  corrugated 
GaAs  crystal  is  shown  in  Fig.  Ill .7.  The  period  is  0.115u  .  Further 
details  on  the  fabrication  are  given  In  Chapter  V. 

In  this  section  we  shall  give  a  short  description  of  the  first 

(5) 

observation  of  a  distributed  feedback  laser  in  GaAs'  .  A  more 
detailed  description  of  the  experiments  is  given  in  Ref.  5. 

The  GaAs  dielectric  waveguide  used  in  the  experiment  was 


■at. 
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produced  by  growing  an  undoped  Ga^AlAs  (x  ~  0.3)  and  a  GaAjs  double 

layer  on  a  GeAs  substrata.  The  thickness  of  the  guiding  GaAs  layer  wan 

2ji  .  The  samples  were  pumped  optically  at  77  *K  using  a  pulsed  dye- 

laser  (Rhodamine  B)  tuned  to  =  63OOA  .  The  individual  pumping 

.0 

pulses  had  a  duration  of  7  x  10  7  sec  and  a  peak  power  of  ~  2hW. 
Cylindrical  lenses  were  used  to  pump  a  rectangular  strip  0.33mm  wide 
and  of  a  variable  length,  as  shown  in  Fig.  III. 8.  The  output  beam  of 
the  laser  emerged  through  a  side  surface  and  was  guided  into  a  spectro¬ 
meter. 

A  spectrum  of  the  laser  output  is  shown  in  Fig.  III. 9.  The 
pumping  intensity  was  1.1  times  the  threshold  intensity  1^  and  the 
length  of  the  pumped  region  was  150p  . 

The  width  of  the  oscillation  spectrum  of  Fig.  II. 9  is  less  than 
lA  and  is  within  the  resolution  of  the  spectrometer.  It  corresponds  to 
a  single  longitudinal  mode  oscillation.  By  increasing  the  pumping 
intensity  and  the  length  of  the  pumped  region,  multi-longitudinal-mode 
oscillation  was  observed  as  shown  in  Fig.  III. 10.  The  length  of  the 
pumped  region  was  700p.  ,  and  the  pumping  intensity  was  1.4  Itfa  . 

The  existence  of  only  one  longitudinal  mode  in  Fig.  III. 9 
rather  than  the  two  symmetrical  lowest  order  longitudinal,  modes  about 
5A  apart,  predicted  by  the  theory  (Fig.  III. 6)  for  a  corrugated 
structure,  may  indicate  that  there  exists  seme  sort  of  periodic  gain 
variation  in  the  sample,  which  contributes  to  the  distributed  feedback. 

In  addition  the  fact  that  close  to  the  center  of  the  spectrum  in  Fig. 

III. 10,  the  modes'  spacings  are  smaller  than  at  the  edges  of  tfr;  spectrum 
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Fig.  III. 8.  Illustration  of  the  laser  crystal  and  the 
configuration  used. 


JJdStx. 


0  "  8240  8250 

Wavelength  (A) 

Fig.  III. 9.  The  oscillation  spectrum  of  a  GaAs  waveguide  laser. 

length  of  the  pumped  region  was  15 Qi  and  i/l^  *= 


supports  the  above  assumption 


As  discussed  in  the  previous  section  the  s pacings  between  the 
nodes  are  quite  similar  to  those  of  a  regular  laser  of  length  L  . 

For  a  regular  laser,  because  of  index  dispersion,  the  longitudinal  node 
spacing  AX  is  given  by 


AX 


2(n  -  X  |£)L 


(III.5.1) 


where  X  is  the  vacuum  oscillation  wavelength,  n  is  the  guide's  index 
of  refraction,  and  L  the  length  of  the  pumped  region.  If  we  take 
n  "  X  fx  =  ^*5  t  a  value  obtained  from  a  Fabry-Ferot  laser  using  a 
similar  waveguide,  we  obtain  AX  ~  lA  .  This  agrees  with  the  spacing 
in  Fig.  1X1.10. 

As  we  have  seen  earlier,  in  a  distributed  feedback  laser  the 
wavelength  of  the  oscillation  is  determined  by  the  period  of  the 
mechanical  corrugation 

X  =  2n2A  (III. 5. 2) 


where  A  is  the  period  of  the  corrugation  and  we  have  assumed 
(t  »  \  /n2)  . 

By  varying  A  it  is  possible  to  tune  the  frequency  of  oscillation. 
A  number  of  waveguides  were  thus  prepared  with  a  different  corrugation 
period.  The  measured  oscillation  wavelength  X  of  the  samples  is 
plotted  as  a  function  of  the  period  in  Fig.  III. 11.  The  tuning  range 
shown  is  about  .  The  mechanical  period  in  this  experiment  is 
three  times  larger  than  the  one  required  by  III. 5.1  and  the  Bragg 


-109- 


reflection  is  provided  by  the  third  order  Fourier  component  of  corru¬ 
gation.  Fig.  111.  11  indicates  that  the  laser  feedback  is  indeed  caused 
by  the  corrugations  and  that  their  use  leads  to  a  stabilization  of  the 
output  wavelength. 


III. 6  The  Merits  of  Distributed  Feedback  Semiconductor  lasers 

The  previous  section  describes  the  operation  of  an  optically 
pumped  GaAs  distributed  feedback  laser.  It  is  evident  now  that  the 
fabrication  of  such  *  ia.<er  is  much  more  difficult  than  that  of  the 
conventional  one  (which  uses  two  cleave  faces  as  mirrors).  Furthermore, 
it  seems  as  if  the  fabrication  of  a  distributed  feedback  injection 
laser  would  be  considerably  more  complicated. 

The  heterostructure  GaAs  lasers,  which  combine  low  threshold 
current  at  room  temperature^  with  relatively  long  life  time^1®^,  are 
multilayer  structures.  The  propagating  mode  is  confined  to  the  inner 
layers  and  does  not  reach  the  surface.  It  is  thus  quite  difficult  to 
perturb  the  mode  by  corrugating  the  surface,  and  another  mean  would 
have  to  be  found.  Corrugating  one  of  the  inner  layers  by  performing 
the  epitaxial  growth  of  the  layers  in  two  separate  steps  is  a  possi¬ 
bility.  However  this  solution  may  tamper  with  the  not  so  well  understood 
life  time  problem  of  these  lasers. 


Because  of  the  above  reasons  we  would  like  to  conclude  this 


chapter  by  discussing  the  advantages  of  the  distributed  feedback  lasers. 
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The  first  advantage  is  the  capability  to  predetermine  the 
wavelength  s' round  which  the  laser  will  oscillate.  As  we  have  discussed 
earlier  the  laser  wavelength,  X  ,  is  given  approximately  by: 

X  «  2  n2  A  (III. 5.1) 

where  A  is  the  corrugation  period  and  n^  is  the  guide's  index  of 
refraction.  It  is  thus  possible  to  vary  the  laser  wavelength  by 
varying  the  corrugation  period. 

The  second  advantage  is  the  discrimination  of  the  distributed 
feedback  mechanism  between  different  longitudinal  modes.  Figure  III. 6 
describes  mode  spectrum  and  required  threshold  gain  for  a  corrugated 
distributed  feedback  laser.  There  are  two  symmetrical  longitudinal 
modes  with  the  lowest  threshold  gain.  The  differences  in  the  required 
gain  between  the  modes  is  a  function  of  the  parameters  of  the  laser 
and  can  be  enhanced.  In  a  laser  where  the  distributed  feedback  is 
supplied  by  gain  variation  tne  discrimination  is  even  better,  because 
there  is  only  one  mode  wi  the  lowest  threshold  gain. 

The  limited  number  of  longitudinal  modes  of  the  laser  is 
important  for  optical  communications  via  optical  fibers.  The  dispersion 
which  exists  in  these  fibers  tends  to  limit  the  information  rate 
capabilities  when  the  bandwidth  of  the  radiation  is  not  narrow  enough. 

In  a  waveguide  laser  the  width  of  its  radiation  depends  in  addition  to 
the  number  of  longitudinal  modes  on  the  number  of  lateral  modes.  Let 
us  therefore,  devote  our  attention  to  these  modes  and  their  control. 
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Tbe  lateral  inodes  of  a  rectangular  waveguide  were  discussed 
in  Chapter  II.  As  an  example  Figure  III.12a  describes  the  (0,0)  (1,0) 
and  (2,1)  modes  of  a  rectangular  laser  (the  mode’s  notation  (m,n) 
is  such  that  mtl  and  o+l  denote  the  number  of  field  extrema  in  the 
x  direction  and  the  y  direction  respectively).  The  thickness  (t) 
of  the  guide  is  usually  l+i  or  less,  while  the  width  (w,1  is  a  few 
tens  of  microns  or  more.  This  means  that  the  phase  constant  in  the 
y  direction  is  much  larger  than  the  one  in  the  x  direction.  In  the 
following  we  shall  see  how  this  large  phase  constant  discriminates 
between  an  (m,0)  and  an  (m,n>0)  modes.  Later  on  we  shall  describe 
a  method  that  discriminates  between  a  (0,0)  and  an  (mX),0)  modes. 
This  will  allow  the  operation  of  a  distributed  feedback  laser  that 
combines  narrow  bandwidth  with  lowest  order  transverse  mode  operation. 

In  a  rectangular  guide  the  propagation  vector  is  given  by  the 
vectorial  sum  of  the  phase  constants  iu  the  three  orthogonal  directions 
x,y,z  .  The  phase  constant  in  the  z  direction  is  p  ,  while  ±  h^  and 
±  h^  are  the  phase  constants  in  the  x  and  y  directions.  The  absolute 
value  of  the  k  vector  in  the  guide  is  thus  given  by 

|k|  =  (p2  +  h2  +  h2)^  (III. 6.1) 

The  vacuum  wavelength  of  the  laser  radiation  \  ,  therefore  is 

2rtn2 

X  ~  (p2  +  h2  +  h2)* 
x  y' 

where  n^  is  the  guide's  index. 


(111.6.2) 


Pig.  III. 12.  (a)  Different  transverse  modes  in  a  rectangular 

cross  section  laser. 

(b)  k  vector  diagram  in  a  distributed  feedback 
laser  of  thickness  t  . 
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As  was  discussed  in  the  previous  sections,  in  a  distributed 
feedback  laser  the  phase  constant  in  the  z  direction  of  the  oscillating 
modes,  is  determined  by  the  corrugation  period.  The  0's  of  the 
lowest  threshold  modes  are  given  approximately  by: 

P  *  i  (111.6.3) 

As  we  have  noted  above  since  t  «  w 


h  »  h 

y  * 


(III. 6. 4) 


and  we  can  ignore  in  III. 6. 2  and  take  m  as  being  0  . 

Using  III. 6. 3,  III. 6. 4  and  III. 6. 2  we  can  express  the  vacuum 
wavelength  as  a  function  of  the  mode  number  (0,n)  . 


(III.6.5) 


(hy)n  is  phase  constant  in  the  y  direction  for  a  particular  mode 
number,  which  for  well  confined  modes  is  given  by 

(hy)n  =  (HI  .6.6) 

Substituting  III. 6. 6  in  III. 6. 5  immediately  yields  the  different 
oscillation  wavelength  of  the  different  thickness  transverse  modes  in 


a  distributed  feedback  laser. 
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(111.6.7) 


As  an  cxaaplc  we  consider  two  transverse  modes  (0,0)  and 
(0,1)  .  The  k  vectors  diagram  for  that  particular  case  Is  shown  in 
Fig.  III.  12b.  The  expressions  for  1(0,0)  and  1(0,1)  cure  given  hy 
III. 6.7  as 


2im„ 


L(o,o) 


far-sr 


L(0,l) 


2iou 

— rm  (III. 6.8) 


Using  typical  values  such  as  A  «  0.115m.  ,  t  =  1m  and  =  3.6  we 
find 

AX 


1(0,0)  ‘  x(o,i)  *  8226  *  80r?QA  c  156A 


This  separation  between  the  two  modes  is  larger  than  half  the 
width  of  the  fluorescence  spectrum  in  Gals!  If  A  is  chosen  properly 
so  that  1^0  corresponds  to  the  center  of  gain  curve,  1(Q 
will  automatically  be  placed  outside  the  gain  curve  and  the  (0,1) 
mode  will  not  oscillate. 

It  is  important  to  realize  that  the  effect  described  above 
is  due  to  the  fact,  that  in  a  distributed  feedback  laser  the  phase 
constant  of  the  longitudinal  mode  is  clamped  to  a  value  of  j  . 

In  a  conventional  laser,  or  the  other  hand,  0  can  take  different 
values  and  thus  allow  oscillation  of  the  (0,1)  mode^11^. 

The  wavelength  separation  between  two  different  width  transverse 
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modes  is  not  large  (because  w  »  t)  ,  and  therefore  does  not  allow 
discrimination  in  the  fashion  described  above.  However  the  wavelength 
separation  is  enough  to  widen  the  bandwidth  of  the  laser  and  should  be 
avoided.  In  addition  the  existence  of  high  order  transverse  modes 
affects  the  far  field  radiation  pattern  of  the  laser  by  causing  the 
beam  to  diverge  strongly.  We  shall  use  another  property  of  the 
distributed  feedback  laser  to  show  how  it  is  possible  to  increase 
the  threshold  gain  required  for  higher  order  width  transverse  modes. 

Let  us  consider  again  equation  III. 2. 11,  which  expresses  the 

value  of  tne  coupling  coefficient  in  terms  of  an  overlap  integral 

between  the  propagating  mode  and  the  perturbation  in  the  dielectric 

constant.  m 

00 

K  e~17|Z  «  £  jj  d£(x,y,z)e^0)(x,y)dxdy  (III.2.11) 

_w 

The  surface  corrugation  thus  perturbs  the  propagating  mode  the  most 
where  it  overlaps  the  peak  of  the  mode's  power  profile.  Figure  III. 13a 
describes  the  power  profile  of  well  confined  (0,0)  ,  (1,0)  and  (2,0) 
modes.  As  can  be  seen,  in  the  (0,0)  case  most  of  the  power  propagates 
close  to  the  center  of  the  laser's  width,  while  in  the  higher  carder 
modes  a  smaller  fraction  of  the  power  propagates  close  to  the  center. 

For  a  high  mode  number  the  power  distribution  tends  to  be  uniform  over 
the  whole  width. 

A  corrugation  extending  only  over  a  central  portion,  say  1/3 
of  the  laser  width,  would  be  felt  mainly  by  the  (0,0)  mode  because 
of  its  unique  power  distribution.  This  will  cause  the  (0,0)  mode  to 


-117- 


have  a  higher  K  and  thus  a  lc^er  threshold  gain  than  the  rest  of  the 
modes.  A  complete  transverse  mode  control  will  thus  be  achieved. 

Let  us  calculate  now  the  degree  of  discrimination  between 
the  modes.  We  first  note  that  the  mode  profile  in  the  x  direction  for 
well  confined  modes  is  given  by 


D,m)<*>  *  (I) 


e,  .(x)  =  (-f  sin 
(0,m)'  '  \w/  w 


x  for  even  m 


x  for  odd  m 


(111.6.9) 


where  w  is  tbo  width  of  the  guide  and  (J^J  is  a  normalization 
factor.  We  can  express  now  the  coupling  coefficient  far  the  different 
modes  in  terms  of  K  and  an  overlap  integral  in  the  x  direction. 

K  is  the  coupling  coefficient  for  the  case  where  the  corrugation 
extends  all  o*ier  the  guide  and  is  given  by  III. 2. 16.  We  can  thus  write 


'  r  1 

K(m,or  KJ  coa2 

-V/2  u 

K(»,of  KJ  8i"2 


for  even  m 


for  odd  m 


(III. 6. 10) 


According  to  III. 6. 10  one  finds  that  for  v  =  j  w 


K(0,0)  -  °-6K  -  *(1,0)  -  °-21C  -  *(2, 0)  ■  °-33K 
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Considering  the  sane  example  discussed  in  section  4  (Fig.  III. 6)  with 

•1 

K  «  6  cm  and  L  »  lmm  ,  we  can  calculate  now  the  threshold  gains  for 
the  lowest  order  longitudinal  mode  hut  with  the  different  lateral  modes. 
The  results  are  shown  in  Fig.  13b  and  indicate  quite  an  effective 
discrimination  between  the  transverse  modes . 

The  last  advantage  of  a  distributed  feedback  laser  which 
should  be  mentioned  is  its  compatibility  with  the  concept  of  integrated 
optical  circuits.  The  conventional  GaAs  laser  is  cleaved  on  both 
ends.  The  high  index  of  refraction  of  GaAs  provided  enough  reflection 
at  the  interface  to  sustain  oscillations.  It  is  difficult,  however, 
to  incorporate  this  method  into  the  fabrication  of  an  optical  circuit. 
One  possibility  is  the  replacement  of  the  cleave  with  a  sharp  step  in 
the  guide's  thickness,  which  will  provide  enough  reflection  for  the 
oscillation.  Such  a  step  was  fabricated  by  ion  milling  and  is  shown 
in  Fig.  III. 14.  The  distributed  feedback  laser  on  the  other  hand  has 
no  discrete  mirrors  at  all  and  therefore  is  highly  compatible  with 
optical  circuits. 

Figure  III. 15a  dt scribes  a  configuration  in  which  a  large 
optical  cavity  GaAs  laser  (LOC)^12^  is  incorporated  into  a  guiding 
structure.  The  guiding  layer  is  the  A1  Ga.  As  wlile  the  active 

Jr  «/ 

layer  is  the  upper  GaAs  layer.  (For  this  purpose  we  must  require  that 
x  >  y  and  that  the  GaAs  layer  be  too  thin  to  support  a  guided  mode 
of  its  own.)  The  gain  is  provided  in  this  case  by  the  optical  pumping 
while  the  feedback  is  supplied  by  the  corrugation.  To  avoid  absorption 
by  the  GaAs  layer  outside  the  laser  area,  it  is  removed  from  there  by 


Bottom 


Fig.  111.14.  A  scanning  electron  micrograph  of  a  3n  deep 
step  ion  milled  in  GaAs. 
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Fig.  III. 15.  (a)  A  configuration  .'or  a  distributed  feedback  large 

optical  cavity  GaAs  laser. 


gradually  tapering  its  thickness  to  zero.  As  opposed  to  the  sharp  step 
discussed  above  and  shorn  in  Fig.  111.14,  in  a  taper  the  mode  gradually 
adapts  itself  to  the  variation  in  the  thickness  and  does  not  suffer 

do} 

severe  reflection  or  loss'  (provided  that  the  final  thickness  guide 
can  support  that  mode). 

Ion  milling,  which  is  described  in  detail  in  Chapter  V,  can 
also  be  used  to  taper  the  width  of  the  laser  and  thus  turn  the  ellipti¬ 
cal  output  beam  into  a  more  round  one.  This  is  shown  in  Fig.  III.  15b. 
The  laser  is  a  distributed  feedback  LOC  laser.  The  corrugation  extends 
over  the  central  portion  only  for  transverse  mode  control.  The  width 
of  the  laser  as  well  as  the  active  layer  are  tapered  to  fit  the  cross 
section  of  a  square  channel  guide.  Such  a  configuration  can  incorporate 
a  single  or  two  longitudinal  modes  together  with  a  round  lowest  order 
transverse  mode  operation. 

It  can  thus  be  concluded  that  the  advantages  of  the  GaAs 
distributed  feedback  laser  encourage  an  attempt  to  implement  it  in  an 
injection  laser  configuration. 
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IV.  FABRICATION  OF  INTEGRATED  OPTICS  COMPONENTS 

IV. 1  Introduction 

All  the  Integrated  optics  components  and  devices  which  we 
have  discussed  above  demand  high  quality  edge  smoothness  and  high 
resolution  formation  in  dimensions  down  to  submicron  sizes.  As  an 
example  we  recall  that  for  a  first  order  distributed  feedback  laser 
in  GaAs  a  corrugation  with  a  period  of  0.115m.  has  to  be  fabricated 
on  the  surface  of  the  GaAs. 

In  collaboration  with  Hughes  Research  laboratories  (Malibu,  CA), 
we  have  employed  ion  beam  milling  and  ion  implantation  in  the  fabrication 
of  the  different  devices.  Following  the  description  of  the  ion  mill Ing 
process  we  shall  detail  a  method  that  was  developed  for  the  fabrication 
of  very  small  period  corrugations  in  a  solid  substrate.  This  method 
utilizes  holographic  exposure  of  photoresist.  The  fabrication  of  channel 
waveguides  (embedded  and  ridged)  as  well  as  directional  couplers  will 
be  discussed. 

Since  in  all  of  these  fabrications  ion  milling  was  used,  we 
devote  the  next  section  to  the  description  of  the  process  and  the 
instrumentation  at  the  Hughes  Research  Laboratories. 
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IV  .2  Ion  Bean  Milling* 

The  sputtering  process  provides  a  convenient  method  for  thin 

film  deposition  and  removal  which  can  be  applied  to  a  vide  variety  of 

materials.  The  term  sputtering  refers  to  the  ejection  of  atoms  from  a 

target  surface  by  boobardment  with  energetic  particles  (usually  ions 

with  energies  in  the  range  of  0.6  to  20  keV).  The  ejected  atoms  can  be 

either  removed  from  a  target  substrate  or  deposited  an  an  opposing 

substrate.  Plasma  sputtering  is  widely  used  in  the  microelectronics 
(2) 

industry'  '  and  has  also  been  used  in  the  fabrication  of  glass  wave- 

(■ i )  i 

guides'  .  The  plasma  is  usually  created  in  an  inert  gas  (at  10  to 

10  Torr)  by  a  dc  or  rf  discharge  and  the  sputtering  occurs  by  ion 

bombardment  from  the  plasma.  Ion  beam  sputtering  is  performed  by 

bombarding  the  target  surface  with  a  collimated  ion  beam  in  high  vacuum. 

This  type  of  sputtering  has  proved  to  be  a  useful  alternative  in  instances 

where  high  vacuum  conditions  are  desired  for  film  deposition,  or  when 

(4) 

control  and  direction  of  micromachining  is  important'  ' . 

For  fabrication  purposes  an  inert  gas  ion  beam  is  produced 
in  a  duoplasmatron-type  ion  source  (as  shown  in  Fig. IV. 1),  accelerated 
through  a  high  vacuum  chamber  and  directed  (with  or  without  final 
focusing)  to  impinge  cm  a  target  surface.  The  high  energy  (typically 
3  to  10  keV)  impact  of  the  ions  causes  sputtering  ejection  of  the 
target  material.  By  means  of  contact  masking  or  by  shadow  masking, 
patterns  can  be  machined  (milled)  into  the  surface  of  the  target. 

Also  referred  to  as  "back  sputtering"  or  "ion  beam  machining." 


Fig. iv. 1  Schematic  diagram  of  the  duopiasmatron 
ion  beam  sputtering  system. 
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The  directed  beam  permits  etching  of  patterns  which  are  more  sharply 
defined  than  those  of  chemical  etching  and,  because  no  material  is 
impervious  to  the  sputtering  attack,  high  resolution  patterns  can  be 
produced  in  a  broad  variety  of  thin  film  materials. 

A  broad  ion  beam  bombarding  a  target  is  used  to  accomplish 
high  vacuum  deposition  of  the  sputtered  target  material  onto  an  opposing 
substrate .  This  is  illustrated  by  the  alternate  substrate  position 
shown  in  Fig. IV. 1.  Electrons  provided  by  the  neutralizing  filament 
are  trapped  in  the  poteitial  of  the  ion  beam  and  are  available  for 
continuous  neutralization  of  any  positive  surface  charge  which  may  tend 
to  accumulate.  Therefore,  dc  bombardment  suffices  for  sputtering 
deposition  or  removal  of  either  metallic  or  insulating  materials. 

Since  the  back-sputtering  rate  proceeds  uniformly  in  a  given 
homogenous  material,  the  three-dimensional  shape  of  a  contact  mask 
will  be  fairly  well  replicated  in  the  substrate  if  sputtering  continues 
until  the  mask  is  just  etched  away.  Some  distortion  of  the  cross- 
sectional  detail  occurs  because  the  sputtering  rate  varies  with  the 
angle  of  incidence.  Also,  the  depth  of  etching  into  the  substrate  may 
differ  from  the  depth  of  the  photoresist  removed  if  the  sputtering 
yields  of  the  two  materials  differ. 


IV.3  Periodic  Corrugation  in  Optical,  Waveguides 


Periodic  structures  in  optical  waveguides  can  be  made  by 


corrugation  of  the  waveguide  surface.  As  we  have  discussed,  when 


choosing  the  right  periodicity  it  is  possible  to  convert  power  from 
one  mode  to  another,  '  or  couple  power  between  a  confined  guided  mode 
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and  an  unconfined  radiation  node,  thus  forming  a  grating  coupler^. 

Mode  conversion  fron  a  forward  propagating  node  to  a  backward 

propagating  node  forma  a  distributed  minor.  This  kind  of  a  mirror 

may  help  In  extending  the  lifetime  of  semiconductor  injection  lasers 

In  the  cases  where  mirror  damage  is  the  dominant  failure  mode.  As  was 

discussed  in  Chapter  III,  if  the  corrugation  extends  along  the  length 

(7 

of  an  amplifying  medium  a  distributed  feedback  laser  results' . 

In  general,  corrugating  the  surface  of  an  optical  thin  film  waveguide 
produces  space  harmonics  with  phase  velocities  which  depend  on  the 
corrugation  period.  These  can  be  used  in  interactions  which  require 
phase  matching  such  as  in  nonlinear  optics  experiments  and  devices 

The  periods  of  the  corrugations  needed  for  these  applications 
vary  from  as  low  as  approximately  O.U41  for  a  distributed  feedback 
laser  in  GaAs  to  tens  of  microns  for  the  phase  matching  of  10. 64m 
second  harmonic  generation  in  GaAs. 

A  new  fabrication  technique  was  developed,  in  which  an 
ion  beam  is  used  to  mill  into  a  substrate  a  grating  mask  produced  by 
laser  holographic  techniques.  The  surface  of  the  guiding  layer  was 
spin-coated  with  a  photosensitive  resist  material  (both  JOPR  or  Shipley 
1350  resists  were  used).  Argon  laser  ( K  =  4580A)  beams  fron  a 
cosnnon  source  were  projected  onto  the  resist  at  an  angle,  as  shown  in 
Fig. IV .2 a.  The  interference  of  these  two  beams  caused  exposure  with  a 
sinusoidal  variation  across  the  surface  at  a  period  of 

^Manufacturers  are,  respectively,  Eastman  Kodak  Company  and  Shipley 
Company,  Inc. 
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By  adjusting  the  exposure  and  development  parameters  the 
pattern  is  left  as  a  partial  coverage  on  the  surface  of  the  substrate 
(Pig.  IV.  2b).  Although  lower  exposure  and  development  leaves  the 
undulating  pattern  on  the  surface  of  the  resist  layer^10^  and,  in 
principle,  the  material  can  be  ion  machined  away  leaving  a  replica  of 
the  undulations  in  the  substrate  surface,  this  approach  is  difficult 
to  control  if  the  resist  material  contains  volatile  components.  This 
is  because  the  ion  beam  heats  the  resist  material  and  gas  evolution 
causes  bubbles  and  blisters  which  lift  the  pattern  and  stretch  it  to 
eventual  destruction.  It  is  more  desirable  to  develop  the  pattern  as 
clearly  defined  and  separated  stripes  which  provides  for  local  gas 
evolution  without  bubble  formation. 

We  have  produced  grating  patterns  with  periods  of  0.28p  , 

0.35n  ,  0.4l+i  ,  1.4+i  ,  and  3. Op  using  the  holographic  technique 

described  above.  Ion  beam  machining  was  then  employed  to  etch  these 
structures  into  the  surface  of  silicon  and  GaAs  substrates.  The  0.4l+i 
pattern  was  etched  to  a  depth  of  0.12+t  into  the  surface  of  GaAs  as 
can  be  seen  from  the  SEM  cross-section  photograph  (Fig.IV.3).  Figures 
IV. 4a  and  b  shew  lower  magnification  photographs  of  the  1.4+i  grating. 
Figure  IV. 5  shows  the  smallest  period  grating  produced  in  the  configuration 
shown  in  Fig.  IV .2 a,  0.28p  .  This  period  was  determined  by  the  laser 
wavelength  rather  than  the  photoresist  resolution.  Equation  IV.3.1 
indicates  that  the  lowest  period  possible  (for  a  =  90  )  is  V2  • 


IV. 3  Cross  section  of  grating  pattern  (0.41  pm  periodicity) 
in  GaAs .  Scanning  electron  microscope  magnification 
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The  smallest  wavelength  continuous  laser  available  is  an  HeOd  (X  =  325C&) . 
Since  even  this  wavelength  is  not  short  enough  for  the  fabrication  of 
the  U0Q&  grating,  the  light  beam  had  to  be  sent  through  a  quartz 
block  (n  m  1.5)  to  further  reduce  the  wavelength.  The  configuration 
used  is  shown  in  Fig. IV .6.  HeCd  laser  beams  incident  on  both  sides  of 
the  quartz  block  reach  the  photo  resist  by  passing  through  an  index 
matching  oil  at  the  bottom  of  the  block.  This  configuration  has 
enabled  us  to  produce  a  115Q&  grating  which  was  shown  in  Fig.  Ill .6 . 
(Special  care  had  to  be  taken  to  remove  the  oil  from  the  photoresist 
before  development  by  dipping  in  xylene . ) 

The  exact  period  of  a  given  grating  was  measured  after  the 
exposure  by  shining  a  laser  beam  on  it  and  calculating  the  period  from 
the  laser  wavelength  and  the  diffraction  angle.  In  the  case  of  the 
115Q&  grating  the  HeCd  beam  must  be  sent  again  through  the  quartz  in 
order  to  obtain  a  diffracted  beam.  Thus  the  index  of  the  quartz 
entered  into  the  calculation  and  it  had  to  be  measured  accurately.  This 
was  done  by  first  measuring  the  period  of  a  grating  with  a  medium  size 
period,  using  HeCd  laser  in  air.  Then  the  third  order  diffraction  from 
this  0.3U65u  grating  was  measured  by  shining  the  HeCd  through  the 
quartz.  From  these  measurements  it  was  possible  to  derive  accurately 
the  index  of  the  quartz  block  (n  =  1.54*0,  and  with  it  the  shorter 
periods  were  determined. 
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Fig. IV. 6  (a)  Holographic  exposure  of  photoresist 

through  a  quartz  block. 

(b)  Ion  machining  of  photoresist  after  development. 


Ion  nachinlng  c cabined  with  either  conventional  photo¬ 
lithography  or  scanning  electron  lithography  is  very  useful  in  the 
fabrication  of  optical  circuitry.  Canponents  such  as  straight  and 
curved  channel  waveguides,  directional  couplers,  thin  film  analogs 
of  lenses  and  prisms etc.  can  be  constructed  by  this  technique. 

The  ion  beam  can  be  used  to  copy  structures  in  the  exposed  photoresist 
(or  electron  resist)  into  the  guiding  material,  or  into  a  specific 
masking  material  that  is  used  whenever  the  guides  are  made  by  ion 
implantation^'  OT  diffusion^ ^ .  We  shall  describe  here  the 
fabrication  of  embedded  and  ridged  channel  guides  in  GaAs. 

Bombardment  of  high  free  carrier  concentration  GaAs  with 
protons  causes  defect  centers  in  the  material,  which  traps  the  free 
carriers  and  thus  produces  a  compensated  layer.  The  thickness  of  the 
layer  is  related  to  the  energy  of  the  protons,  and  for  the  300  keV 
protons  used  it  is  3 u  .  The  presence  of  charge  carriers  leads  to  a 
depression  of  the  optical  dielectric  constant.  The  difference  Ng  -  Ng 
in  carrier  density  between  the  guide  and  substrate  causes  the  dielectric 
constant  of  the  guide  to  exceed  that  of  the  substrate  by: 

(Na  -  N )e2 

Ae  -  —  •  »  | -  (3V.4.1) 

m  cu 

where  N  and  N  are  the  carrier  densities  in  the  substrate  and  in 
s  g 

# 

the  guide  respectively,  m  and  e  are  the  effective  mass  and  charge 
of  the  carriers  and  oo  is  the  radian  optical  frequency.  As  we  shall 
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see  later  this  plasm  contribution  seeas  to  dominate  for  wavelength 
away  from  the  bandgap  of  the  material.  Since  a  gold  layer  1.5m  thick 
can  completely  block  the  300  keV  protons  used,  «n  optical  circuit  can 
be  implanted  if  a  suitable  gold  mask  is  constructed  on  the  surface  of 
the  GflAs  substrate. 

In  the  fabrication  of  channel  waveguides  closely  spaced  so 
as  to  allow  directional  coupling  between  them,  we  have  used  a 
conventional  photographic  mask.  A  Gafts  substrate  was  sputter  coated 
with  a  1.8m  thick  gold  layer.  Finally,  a  photoresist  layer  (2m) 
was  spin  deposited  and  was  exposed  through  a  mask  of  appropriate 
transparent  stripes.  The  exposed  resist  was  developed  down  to  the  gold, 
to  form  stripes  of  resist  with  a  ridge  shape.  The  succeeding  step  was 
ion  machining  of  the  clear  areas  of  gold  down  to  the  substrate.  Figures 
II. 4a  and  II. 4b  show  the  resulting  gold  stripes  with  the  remainder  of 
the  photoresist  on  top  of  them.  The  distance  between  the  centers  of 
adjacent  channels  is  about  6.3m  •  The  clear  areas  between  the  gold 
stripes  are  where  the  proton  implantation  takes  place  and  the  channel 
waveguides  are  formed.  After  the  implantation  the  gold  was  removed  and 

O 

the  samples  annealed  at  500  C  for  10  minutes  in  order  to  reduce  the 
optical  attenuation.  The  optical  measurements  of  these  guides  and  the 
coupling  between  them  is  described  in  Chapter  II. 

One  of  the  important  parameters  that  influences  these  measure¬ 
ments  is  the  dielectric  discontinuity  between  the  substrate  and  the 
Implanted  region.  An  estimate  for  this  dielectric  discontinuity  is  the 
plasma  contribution  given  by  (IV.4.1).  However  according  to  it,  the 
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change  in  the  index  of  refraction  at  X  =  1.06u  should  be  smaller  than 
at  X  s  1.15u  .  The  measured  values  of  the  coupling  coefficient  K 
(see  Chapter  II)  contradict  that,  and  indicate  an  increase  in  Ae 
at  1.06fi  rather  than  a  decrease  with  respect  to  1.15u  .  Because  of 
that,  an  attempt  was  made  to  measure  and  compare  the  values  of  the 
dielectric  discontinuities  at  both  wavelengths.  Since  it  is  quite 
difficult  to  measure  small  changes  in  the  index  of  refraction  that  take 
place  in  such  small  dimensions  the  following  method  was  used.  The 
multichannel  sample  shown  in  Fig.  II. 4b  was  polished  from  the  substrate 
side,  and  used  as  a  phase  grating  by  shining  the  laser  light  perpendicular 
to  the  channels .  The  geometry  and  the  resulting  diffraction  orders 
are  shown  in  Fig.  1V.7.  Using  the  amplitudes  of  the  diffraction  orders 
for  1.15u  and  1.06u  it  is  possible  to  find  the  ratio  between  An^  ^ 
and  An^  ^  .  It  was  found  that 

An. 

"  ■■■'?=  «  1.1  (iv.4.2) 

An.  .c 

1.15u 

This  means  that  as  the  wavelength  approaches  the  bandgap  of  GaAs  (0.9n) 
the  plasma  contribution  is  not  the  dominant  contribution  to  An  , 
and  that  additional  contribution  (probably  because  of  increased  bandgap 
tail  absorption  in  the  implanted  region)  starts  to  increase  An  . 

In  order  to  derive  the  absolute  change  in  the  index  from  this  measure¬ 
ment  an  estimate  of  the  channel  dimensions  has  to  be  made.  Using  the 
mask  shape  as  shown  in  Fig.  11.4b  (width  separation  3.9-0  and 
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a  penetration  depth  of  34  for  the  protons,  the  approx inn te  values 
for  An  are: 

4n1.06  *  5-5'10'3  4”l.l5  ■  5'10'3 

while  the  value  calculated  using  the  plasma  contribution  IV. 4.1  for  a 

1 A  o 

substrate  concentration  of  N  =  2.10  cm"-5  are 

8 

4n1.06  "  3-7'10'3  '  4-3-10'3 

The  agreement  thus  is  adequate  for  1.154  only. 

As  can  be  seen  in  Fig.  II. Ub  the  side  walls  of  the  gold 
stripes  are  quite  rough.  The  reason  is  the  limited  quality  of  the 
photographic  mask.  A  better  wall  definition  can  be  achieved  by  using 
one  of  the  methods  described  earlier,  namely,  scanning  electron  litho¬ 
graphy,  or  the  holographic  technique.  The  fact  that  the  attenuation 
of  the  channel  guides  is  about  that  of  planar  guides  fabricated  in  the 
same  substrate  (about  5  cm"1)  indicates  that  the  walls  are  smoothed 
out  by  the  nature  of  the  proton  implantation  and  the  subsequent 
annealing .  However,  rough  walls  are  a  major  limitation  in  the  fabrica¬ 
tion  of  ridged  guides.  Therefore  in  the  fabrication  of  these  guides 
(which  is  shewn  in  Fig.  IV .8)  a  mask  produced  by  the  holographic 
technique  was  used  to  expose  the  photoresist  on  top  of  the  GaAs. 

Towards  the  end  of  the  machining  the  sample  was  tilted  at  various  angles 
with  respect  to  the  ion  beam  to  ensure  the  effect  of  "ion  polishing." 
Two  examples  of  ridged  guides  are  shown  in  Fig.  11.1$.  The  rough 
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Fig,  IV.8  Channel  waveguide  fabrication  by  removal  of 
superfluous  sections  of  an  epilayer  by 
ion  filling. 
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channel  was  fabricated  with  a  conventional  nuk  and,  as  can  be  seen 
in  the  photograph,  with  a  photoresist  which  was  not  thick  enough  to 
protect  the  top  of  the  channel.  Figure  II. 15b  on  the  other  hand 
describes  a  guide  fabricated  with  the  use  of  a  holographic  mask  and 
an  "ion  polish"  at  the  end  of  the  machining.  The  smoothness  of  the 
walls  in  this  figure  as  well  as  of  the  step  in  Fig.  Ill .  14  indicates 
that  ion  milling  can  be  used  to  fabricate  optical  resonators  for 
surface  lasers  and  terminations  for  waveguides  (in  materials  that  do 
not  cleave  too  well)  for  input,  output  and  coupling. 


-143- 


Refercnces  for  Chapter  IV 

1.  H.  L.  Garvin,  E.  Garmire,  S.  Sanekh,  H.  Stoll,  and  A.  Yariv, 

"Ion  Bean  Micronach lning  of  Integrated  Opt.’cs  Components, "  Appl. 

Orptica  12,  455  (1973). 

2.  Many  applications  have  been  cited  at  the  Symposia  on  the  Deposition 
of  Thin  Films  by  Spattering,  sponsored  by  the  Univers  ity  of 
Rochester  and  Bendix  Corporation,  Rochester,  New  York,  held  in 

19 66,  1967,  and  1969,  and  at  the  Conference  and  School  Sessions  on 
"The  Elements,  Techniques  and  Applications  of  Sputtering, "  sponsored 
by  the  Materials  Research  Corp.,  Orangeburg,  New  York. 

3.  J.  E.  Groell  and  R.  D.  Standley,  "Sputtered  Glass  Waveguides  for 
integrated  Optical  Circuits, "  BSTJ  48,  3445  (1969) . 

4.  H.  L.  Garvin,  "Ion  Beam  Sputtering  for  Thin  Film  Deposition  and 
High-Precision  Mlcrcmachlning, "  Proc.  3rd  Symposium  on  the  Deposition 
of  Thin  Films  by  Sputtering,  Rochester,  N.  Y.,  pp.  4-11,  September 
7-10,  1969. 

5.  D.  Marcuse,  "Mode  Conversion  Caused  by  Surface  Imperfections  of  a 
Dielectric  Slab  Waveguide,"  BSTJ  48,  3177  (1969). 

6.  M.  L.  Dakss,  L.  Kuhn,  P.  F.  Heidrich,  and  B.  A.  Scott,  "A  Grating 
Coupler  for  Efficient  Excitation  of  Optical  Guided  Waves  in  Thin 
Films,"  Appl.  Phys.  Lett.  16,  523  (1970). 

7.  H.  Kogelnik  and  C.  V.  Shank,  "Stimulated  Emission  in  a  Periodic 
Structure, "  Appl.  Phys.  Lett.  18,  152  (I97I). 

8.  S.  Sanekh  and  A.  Yariv,  "Phase -Matchable  Nonlinear  Optical 
Interactions  in  Corrugated  Thin  Films,"  Appl.  Phys.  Lett.  21,  140  (1972). 


IMUH 


-144- 

9.  H.  W.  Yen,  M.  Nakamura,  E.  Garmire,  S.  Somekh,  A.  Yariv,  and 

H.  L.  Garvin,  "Optically  Pumped  GaAs  Waveguide  Lasers  with  a 

Fundamental  0.12|i  Corrugation  Feedback, "  to  be  published  in 
Optics  Communications. 

10.  R.  G.  Braudes  and  R.  K.  Curran,  r  Modulation  Transfer  Function  of 
AZ111  Photoresist,"  Appl.  Optics  10,  2101  (1971). 

11.  R.  Shubert  and  J.  H.  Harris,  "Optical  Guided-Wave  Focusing  and 
Diffraction,"  JGBA  61,  154  (1971).  Also,  R.  Ulrich  and  R.  J. 
Martin,  "Geometrical  Optics  in  Thin  Film  light  Guides,"  Appl. 
Optics  10,  2077  (1971). 

12.  E.  Garmire,  H.  Stoll,  A.  Yariv,  and  R.  G.  Hunsperger,  "Optical 
Waveguidlng  in  Proton- Implanted  GaAs,"  Appl.  Phys.  Lett.  21, 

87  (1972). 

13.  R.  D.  Standley,  W.  M.  Gibson,  and  J.  W.  Rodgers,  "Properties  of 
Ion-Bombarded  Fused  Quartz  for  Integrated  Optics,"  Appl.  Optics 

u,  1313  (1972). 

14.  W.  E.  Martin  and  D.  B.  Kali,  "Optical  Waveguides  by  Diffusion  in 
II-VI  Compounds,"  Appl.  Phys.  Lett.  21,  325  (1972). 


Mi 


ft'-. 


-145- 


Appcndlx  I;  The  Modes  of  a  Planar  Waveguide 

The  modes  of  the  waveguide  are  derived  by  solving  the  wave 
equation  and  requiring  the  transverse  components  of  the  electric  and 
magnetic  field  to  be  continuous  across  the  waveguide  discontinuities. 

In  the  case  of  the  planar  guide  there  are  two  kinds  of  solutions,  IE  in 
which  the  electric  field  is  transverse  and  1M  in  which  the  magnetic 
field  is  transverse.  For  a  given  frequency  the  guide's  thickness  and 
dielectric  discontinuity  determines  the  number  of  guided  modes.  These 
guided  modes  peak  inside  the  guide  and  decay  exponentially  outside  the 
waveguide  region.  A  different  sort  are  the  so-called  "radiation  modes” 
which  do  not  decay  outside  the  guide,  but  extend  from  x  =  to 
x  *  +®  and  are  capable  of  carrying  power  away  from  the  waveguide^ . 

Let  us  describe  now  the  different  guided  modes  of  the  geometry  shown  in 
Fig.  A. la. 


TE  Modes 


The  field  components  of  these  modes  are  Ey,  and  The 


field  component  E  obeys  the  wave  equation 


2  ni2  *2Ev 
^  E  *  — — — i  i  =  1,2,3 

y  at2 


(A.I.l) 


We  take  £  (xzt)  in  the  form 

y 


E  (xzt)  =e  (x)  e1(a)t"pz) 

y 


(A. I. 2) 


«i2  - •»  propagation 


(b) 

Fl«.  A.l  (a)  The  basic  configuration  of  a  planar 
dielectric  waveguide. 

(b)  Mode  dispersion  in  the  above  waveguide. 
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*  I W  *  &sr_J  I  Vx)  *2dx  ■ 1 


(A.I.7) 


Using  A.  1.3  ve  determine 


C  2h  — 


P(tf|  ♦  i)(hc+q‘) 


(A. 1.8) 


1M  Modes 


The  field  components  of  these  modes  are  H  ,  E  and  E 

y*  x  z 

H  (xzg)  -  *(x) 

J 


^(«t)  .14.1  „(x)ei(<ot-«>.) 
<»  Sz  at 


^(xzt)  =  -  —  ^ 
«*  dx 


The  mode  profile  M(x)  is  taken  as 


C  S  e”** 


c|jr  cos(hx)  -  sin(hx)j 
cjr  cos(ht)  +  sin(ht)J 


0  *  x  < 


■t  S  X  s  0 


-®  <  X  s  -t 


(A.1.9) 


The  eigenvalue  equation  is 


ht  •  (mfl)«  -  tan"1  -  -  tan”1  $ 

P  <1 


(A.I.10) 


IWLa  condition  determines  the  value  of  C  as^2^ 


From  the  above  solutions  it  is  evident  that  the  propagation 
constant  in  the  z  direction  -p  of  a  guided  mode  varies  between  the 
two  limits 

n3k0  *  P  S  n2k0  (A. I. 12) 

(We  have  assumed  that  n3  >  ^  which  is  usually  the  case.) 
p  is  less  than  n2kQ  in  order  to  allow  the  cosine  or  sine  variation 
along  the  x  direction  inside  the  waveguide.  Outside  the  waveguide  the 
mode  profile  should  fall  exponentially  and  this,  according  to  A. 1.4,  means 
P  >n3kQ  .  The  exact  value  of  p  as  we  have  seen,  is  determined  for  a 
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given  frequency  by  the  guide's  thickness  and  the  dielectric  discontinuity. 
At  cutoff  i>  equals  n^k^  and  as  the  frequency  inc^e^ses  0  approaches 
n^kQ  .  This  behavior  is  illustrated  in  Fig.  A. lb,  for  propagation  in 
the  ±z  directions . 

Also  described  in  this  Figure  A. lb  is  a  continuum  of 

radiation  modes  whose  0  is  less  than  n^kQ  .  This  means  that  the  p 

P  2  2  h 

of  these  modes  is  Imaginary  (p  *  (O'  -  n^  kQ  )*)  and  therefore  the 
node  does  not  decay  exponentially  in  the  x  direction  but  rather 
oscillates  from  to  +“  .  Any  mechanism  that  couples  power  from  a 
guided  mode  into  a  radiation  mode  will  thus  cause  a  power  leak  from 
the  guide^. 
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Appendlx  II:  The  Relative  Signs  of  the  Coupling  Coefficient 

Let  us  consider  two  lossless  modes  a  and  b  ,  whose 
amplitudes  A  and  B  are  related  by: 

i  -  -V +  V  (*•»•« 

i  -  -«bB  +  VA 

We  consider  first  the  case  where  both  modes  carry  power  in  the  same 
direction.  The  conservation  of  total  power  at  any  given  z  is 
expressed  by 

~  =  (AA*  +  HB*)  =  0  (A.II.3) 

where 

Pa  =  AA*  Pb  =  BB*  (A.II.4) 

Because  of  the  plus  sign  in  A.II.3,  it  is  satisfied  when 

-  -C  (A-M-5) 

Next  we  consider  a  case  in  which  the  two  modes  carry  power  in 
opposite  directions.  The  conservation  of  power  is  described  now  by: 

—  (AA*  -  SB*)  «  0  (A. II. 6) 

Because  of  the  minus  sign  the  equation  is  satisfied  only  when: 

-C 


(A.II.T) 


